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This article is devoted to the asymptotic expansion of the generalized 
Harish Chandra-Itzykson-Zuber matrix integral for non-unitary symmetries 
characterized by a parameter (3 (as usual /3 = 1, 2 and 4 correspond to the 
orthogonal, unitary and symplectic group integrals). 

The results are of the form J2perm. exp(X) XiK)f{'Tij)/[Y[i<j ''"ij]^"^; iii which 
Xi and Aj are the eigenvalues of the two kxk matrices, and Tij = {xi — Xj){Xi — 
^ ■ Xj) . A WKB-expansion for / is derived from the heat kernel differential 

equation, for general values of k and (3. From an expansion in terms of zonal 
polynomials, one obtains an expansion in powers of the r's for (3 = 1, and 
generalizations are considered for general f3. A duality relation, and a trans- 
formation of products of pairs of symmetric functions into r polynomials, is 
used to obtain the expression for /(tjj) for general j3. 
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1 Introduction 

We study the HIZ-integral / 

/(X,A) = /"exp(tr gXg'^ k)dg 
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where g is one of the Lie groups 0{k), U{k) or Sp{k), and dg is the corre- 
sponding Haar measure. 

The k X k matrices X and A are Hermitian (real, complex and quater- 
nion). This integral appears at several key points in random matrix theory. 
When g varies over the unitary group, the result is well-known : due to 
Harish- Chandra [11], it has been rederived in the context of random matrix 
applications by Itzykson and Zuber [12]. The HIZ result turns out to be 
"WKB exact" in the unitary case, i.e. to be equal to the sum of Gaussian 
integrals over the k\ saddle-points of the integrand. 

For non-unitary cases, this semi-classical property is no longer true. How- 
ever the full WKB expansion is needed in several problems arising in random 
matrix theory, for instance in the study of universal short-distance correla- 
tions of the eigenvalues. There are several studies of those non-unitary cases. 
For instance, in the orthogonal case the result is known as an expansion 
in terms of zonal polynomials Zp{xi, ...,Xk) [15], which are symmetric poly- 
nomials in the Xi. However, the expressions in terms of zonal polynomials 
do not provide the desired WKB expansion around the saddle-points, which 
involves combinations of the form Y\i^j{xi — Xj). Furthermore there is no 
explicit expression for zonal polynomials of arbitrary order. 

For general f3, zonal polynomials are generalized, and known under the 
name of Jack polynomials [13] but the situation is not improved. In this 
article we extend the HIZ result to non-unitary groups and obtain expressions 
of the form 



e^^'''x(r.,)/n 






or 

where Tjj = {xi — Xj){Xi — Xj), with the parameter (3 corresponding to the 
0{k),U(k) and Sp{k) group integrals for /? = 1,2 and 4, respectively. Al- 
though the group integrals correspond to those three values, the expressions 
may be continued to arbitrary j3 through the differential equation for /(X, A). 
The fact that, once the exponential prefactor e^^'^* is extracted, the result 
may be expressed in terms of the variables r's is far from trivial. Although it 
follows from this paper, through a first representation in terms of extended 
zonal polynomials, and further identities on products of symmetric functions 
in the variables Xi and Aj, an a priori proof would be welcome. 



Recently, a similar expansion, using the formalism of Baker-Akhiezer 
functions, has been considered for /(X, A) by Berest [2]. 

However this is still not in terms of the Tij, and it involves differential 
operators which make the calculation of /(X, A) difficult, except in simple 
cases. 

In the case of the symplectic group (/3 = 4), we have found earlier that the 
WKB expansion for /(X, A) terminates after a finite number of terms[4, 5]. 
We have obtained these explicit finite expansions for k=2,3 and 4. In this 
article we study the coefficients of this expansion for general k ; however the 
structure for arbitrary k is still only partially known. Again a proof of the fact 
that the WKB expansion terminates after a finite number of terms would be 
welcome. Recently the problem has been considered by Ben Said and 0rsted 
[1] but they have simply verified like ourselves this amazing property. 

For other even integer values of (3, {(3 = 6, 8, 10, ...), similarly the WKB 
series stops after a finite number of terms. In this article we discuss also the 
coefficients of these expansions, but again proofs of the finiteness are lacking. 

The expressions that we have obtained are of the form 

where / is given as a series in terms of the Tij = (xj — Xj){\i — \j), 

f(r,j) = l+c,J2T,, + 0ir'). (3) 

i<j 

and the successive coefficients Ci,C2,--- are functions of [3 and k. The 
"perm." in (2) means a sum over the k\ permutations of the Aj's. 
In the case k = 3, P = 4, we have 

f{Tij) = 1 - -(ri2 + T23 + Tig) + --{T12T23 + T23T13 + Ti^Tu) 

- — T-12T-23n3- (4) 

In the case k = 4, (3 = 4, the result is again a similar finite polynomial which 
is given in [4]. Our aim is to find similar compact expressions for general 
values of k and /3. 

Let us note that in the cases f3 = 2m{m = 2,3,...), since the function 
f{Tij) is a polynomial, it means that the HIZ integrals may be written as 

^/3 = E e^'""' ^ f[1 + «i E — + 0(A)] (5) 

pi^. [A(a:)A(A)]f ^ n^ \^" 

which is a WKB expansion corrected by a finite number of terms. 



The integral / is also known to be expressible as an infinite sum over 
extended zonal polynomials Zp{x), with a parameter « = f [15, 13, 16] 



'/3 






m=0 



where Xp(1) is a character, and Zp{T) is a dimensional constant which depends 
upon k. (Explicit values of Zp{x), Xp(l)) Zp(l) are given in Appendix A for 
the lower orders). 

To extract the WKB-exponential e^^'^' factor in (2), we shift X ^ X = 
X — -^trX. Then, the above expression becomes 

J E->A.-iE.<,-. Y 1. I y-^ Zpix)ZpiA) 

In the case f3 = 1, the power of the Vandermonde factor in the denomi- 
nator of (2) vanishes. Therefore one may compare directly the expression of 
(7) with (2) ; consequently the expression for /(tjj) is obtained by rewriting 
the products of symmetric functions in terms of polynomials in the r's. 

In the other cases, (3^1, one needs to extract the Vandermonde fac- 
tor. For instance if f3 is an even integer, this factor is antisymmetric under 
permutations of the Xj and it makes the problem difficult. 

In this article a dual representation is derived , which solves the problem 
of the Vandermonde factor. The duality comes from another expression of 
the integral Ip in terms of extended zonal polynomials. Namely, we can write 
the series (7) as 

perm. i.ii<j[{^i ~~ ^j){^i ~^ ^j)\ 
_iV r- ^ 1 1 ^ . ,ZJX)ZJA), 

where the parameter a is now ^^ instead of |. This is a duality transforma- 
tion of /9 — > 2 — /3. Note that it is the dual invariant product /3(| — 1), which 
appears as a coefficient of the inverse square potential in the Calogero-Moser 
model. The expression (8) is close to the desired expression (2). The WKB- 
exponential term e^^'^* and the correct power of the Vandermonde term are 
already present in (8). The two expressions (7) and (8), are dual under the 
transformation (3 ^ 2 — (3. Then if one writes the subsequent series in terms 
of the variables r, one obtains the desired WKB expansion. 
This article is organized as follows: 



In section two, we discuss the differential equation satified by the integral 
/(X,A). 

In section three, the series expansion in terms of Tj^ = (xj — Xj){\i — \j) 
is obtained by consideration of residues (residual equations). 

In section four, the case /9 = 4 is investigated separately; the r-expansion 
of / is discussed. (The result up to fourth order is given in Table A). 

In section five, we extend the r expansion to arbitrary j3. The residual 
equations are also used at this effect. (The coefficients C[graph] are given in 
Table B). We find that the residual equations have no P dependence. 

In section six, we derive the r expansion directly for [3=1 from the zonal 
polynomial expansion of the integral /, based on the expression (7). ( These 
results are expressed in Table C). 

In section seven, we extend the result of the section six to general (3, 
using the extended zonal polynomials (Jack polynomial). We find a duality 
relation of the type (8). By the duality trasformation a = | — ;> ^^, the 
expressions for C[graph] given in section five (Table B) are rederived. The 
residual equation of section five, is shown to be also valid for the coefficients 
of C[graph], which is determined from the extended zonal polynomials. 

In section eight, we discuss the case /3 = 4 by taking a = —1 in (8). In 
the denominator of (8), one has to deal with the divergent factor j^. This 
divergence is cancelled by the sum over the partition of the zonal polynomials. 
We explicitly evaluate (8) up to order six in the limit a — > — 1. We find the 
large k behavior of this expansion. 

In section nine, the case (3 = 2, where the duality transformation becomes 
singular, is discussed. 

In section ten, the large k limit and large f3 limits are discussed briefly. 

Finally, a summary and discussions are given. The application of this Tij 
expansion to the random matrix theory is briefly mentioned [10]. 

In Appendix A, the extended zonal polynomials in terms of the classical 
symmetric functions s„, a sum of powers, are given up to order six for general 
a. The characters Xp(l) and the dimensional constants Zp{T) are explicitly 
given. 

In Appendix B, the expression of paired products of classical symmetric 
functions, s„(x)s„(A), as polynomials in r , is investigated with the help of 
differential operators Dj^j^y,','., which act on r series. This technique is used 
to prove various identities which appear in the r transformations. 

In Appendix C, we discuss some cubic and quartic identities for the vari- 
ables Tij. In view of these identities, the r expansion has some freedom which 
is used to fix some of the coefficients in the expansion arbitrarily. 

In Appendix D, we characterize each r terms by the unique factors 
x'llx'll ■ ■ ■ AJJA^2 ■ ■ ■) which are picked up by the differentials D*^*^'" . Using 
these unique characteristic differentials, we obtain the coefficients of r ex- 



pansion for /. 

2 Differential equations for the HIZ-integrals 

The Laplacian operator with respect to the matrix elements 

L=i:^- (9) 

(a short-hand notation for the appropriate group invariant Laplacian) acts 
on the integrand of I{x, A) in (1) by 

^gtrAX _ (trA2)etrAX_ (^g) 

Given that I{x, A) is a function of the eigenvalues Xi of the matrix X, the 
equation (10) reads [3], 

fc Q2 k Id 



with the eigenvalue e 



e = j:\l (12) 



Note that the integral / is manifestly symmetric under interchange of the 
matrices A and X, but the procedure is dissymetric. The solutions will of 
course restore this property, which is far from obvious if one considers the 
equation (11) alone. The x-dependent eigenfunctions of this Schrodinger-like 
operator have a scalar product given by the measure 

{(pi\(p2) = / dxi ■ ■ ■ dxk\A{xi ■ ■ ■ Xk)f (flixi ■■■Xk)ip2{xi---Xk) (13) 

The measure becomes trivial if one multiplies the wave function by |A|^/^ . 
Thus if, for some given ordering of the Xj's, one changes I{x) to 

ilj{xi---Xk) = \A{xi---Xk)f/^I{xi---Xk), (14) 

one obtains the quantum Hamiltonian, 

This Schrodinger equation is a simple Calogero-Moser model. 



For a given ordering of the Xj's, pulling out of (15) the exponential factor 
which is the value of the integrand at its saddle-point, one obtains 

^Pq = e^i^i+'-'+^fe^fe;^ (16) 

where x satisfies 

[E^ + 2EA.^-yE7^47^]x = (17) 

in which y stands for 

l/ = /5(f-l). (18) 

In [5] an expansion of x for large r's 

x = i-|(E-^) + |(-i + |)(Ei) + --- (19) 

has been introduced. 

For even integer values of [3, this expansion terminates after a finite num- 

k 

ber of terms. The term of highest degree is [JJ'Tij]'^ . Therefore, one 

i<j 

may write the whole expression as 

From the definitions of ip and x ^^ (14,16), the integral / is then expressed 
as 

^/^(^'^)= E rn l^,-i ^^-'^'''f^ (21) 

perm.ofAi Ulj<i '«iJ 

in which we have normalized fp such that /^ = 1 + 0{Tij). 

Continuing from even integer /9's to arbitrary real numbers, we consider 
below the expression of //? for arbitrary [3 and arbitrary k. 

3 Expansion in powers of r 

We now perform an expansion in powers af the variables r. Let us begin with 
the simple case (3 = 4, and k = 3 . The differential equation for ip, defined 
by (15), is 

P^ = (22) 



where the operator P is 

P = J2-f^-J2jzr^zr^-J2^l (23) 

The expression of ip in terms of Tij consists, in this simple k=3 case, of the 
sum of four terms of increasing degrees: 

^ = V^O + ^1 + ^2 + ^3 (24) 

^2 = Cae^^'^'-— — — — [ri2r23 + T23T31 + r3iri2] (27) 

A{x)A[A) 

^3 = ^3eS-^'^^^[r,2r23r,3] 

= Cge^^^^^ (28) 

where A(a;) = {xi - a;2)(a;i - a;3)(a;2 - X3) and ri2 = {xi - a;2)(Ai - A2). 
The unknown coefficients C„ may be obtained from the differential equation. 
Applying P on ipo and ipi, we have 

A(a;)A(Aj Xi — a;2 Xi — x^ X2 — X3, 



p^^ = C'ieS^^^'[- 



(^3 - a;i)(x2 - a;3)ri3r23 {xi - 3:2) (x2 - a;3)ri2r23 
2 4 4 , 



{xi - a;2)(a;3 - a;i)ri2ri3 {xi - 0:2)^13X23 (a;i - X3)2ri2r23 
4 2 2 2 



(a;2 - a;3)2ri2ri3 (xi - 0:3)2x23 (a;2 - X3)Vi3 {xi - ^2)^23 

2 2 2 

{x2 - a;3)2ri2 (xi - ^2)2^3 {xi - xs^tu 



where the ffrst three terms are generated by -p^, and the next three ones 



i! 



from the factor —4 T^ in P. The last six terms are obtained by 



the differentiations both of exp[X]Aja;j] and — . Note that Pipi consists of 

Tij 

terms of order 0(l/a;^) and 0{l/x^). 

The unknown coefficients C„ are determined by the condition Pip = 0. 
Since Xi and Aj are arbitrary variables, we have the freedom of taking any 
particular limit ; for instance xi -^ X2- From the condition Pip = 0, 
the pole term of l/(xi — X2) present in Pipo should be cancelled by the 
term coming from Pipi. This requirement is sufficient to fix the coeffi- 
cient Ci. Indeed In Pipo, the residue of this pole is proportional to Ai, 

Pihr, ~ —2- , , , X » /^ X ■ From the second and the fourth terms in 

(a;i - a;2)A(a;)A(A) 

Pipi, we have — 6C1- . . , . . ., . term. Therefore, the required can- 

[Xi - X2)/\{X)A{\) 

cellation fixes the coefficient Ci to be —1/3. 

C - 

The term of —2- — in Pipi is cancelled by the term in P'?/'2- 

[Xi - X2YT23 

Indeed we have 



p^2 = C2eS"»^'[- — ( ^ ^ 



ri2 {xi - x^Y {X2 - X-i\ 



2' 



4,1 1,4,1 1 

+ 7 T77 7 777 + 



T-23 {Xi - X^Y (yXi - X2Y ri3 (Xi - X2y {X2 - X^f 

2 2 2 , 



(Xi - X2y {Xi - X^Y {X2 - X^] 



2J 



(31) 



C 
The pole term —2- — should be cancelled, and thus 

(Xi - X2)^r23 

2Ci + 4C2 = 0. (32) 
Thus one obtains C2 = 1/6 . The last term ■?/'3 is 

V'3 = C73eS"''' (33) 

Pi', = Qei: -'"'(- ^ - , ^ ., - , ^ ,, ) (34) 

(Xi - X2)^ (X2 - X3)^ (Xi - X3)^ 

The cancellation of the double pole l/(xi — X2)^ requires 

2C2 + 4C3 = (35) 

i.e. C3 = -1/12. 

We have thus determined the coefficients of the various terms in the 
expansion for /5 = 4 , /c = 3. The results obtained in this manner coincide 
with our earlier results [4] ; this justifies the expression (4) given in the 
introduction. 
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The remarkable point is that the cancellation conditions connect linearly 
the term C„_i and the n-th order term C„ as 

ACn-l + BCn = 0. (36) 

It turns out that this structure holds for arbitary k and f3, and it turns out 
that the coefficients A and B hereabove are functions of k and not of (3. 



4 The r-expansion for the symplectic case (3 = 

4. 

The group parameter j3 remains equal to 4 in this section. (In the next 
section, we will consider general (3). Generalizing the simple k = 3 example 
of the previous section, we examine the cancellation of the pole l/(a;i — X2) 
when Xi -^ X2 for arbitrary k. Let us focus on a specific pole term in the 
various pieces of Pip. 

(i) In Pipn-i, the action of 7—5- + 7—^ on e^'^'^^ and l/(a;i — X2) yields 

oxf 0x2 

x^-X2^ A(x)A(A) -^^ '™^- 

(ii) In Pipn, the multiplication by —4- — gives 

(xi -X2y 

xi-X2^ A(a;)A(A) -l^ '™^- 

£ 1 
(iii)The derivative — ^ of ipn yields —2- (where m 7^ 

dx2 [Xi - a;2)(a;2 - Xm) 

1,2). This term may be written as 

g^A.:r. r_2- ^ 7 n^u]- Si'^^^ Tim/ix2 - x„,) rcduces 

A(a;)A(A) {x^ - X2){x2 - Xm) 
to Ai — Am in the limit x\ —^X2, we are left with the pole term 

_^ Ai-Am ^^A,. ^\^^ nr„- in this limit. 
xi-xs A(a;)A(A) ^^ ^ 

Va 1 

(iv)Leaving aside the overall factor e^ '^'-t^^-tttti the pole terms in 
^ ^ A(a;)A(A) 

(i),(ii) and (iii) must cancel. The residue of l/{xi — X2) is the combination 

of Ai n Tij with unknown coefficients C for each monomial in ipn- 

For instance, in the case of k=3, n=l, from (i),(ii) and (iii), we would 

find 

-2Co^^ - 4Ci^^ - 2Ci^^ = (37) 

Xi — X2 X\ — X2 X\ — X2 

and the cancellation imposes Co + 3Ci = (Co = !)• The pole term in 
Pip is always a single pole, once we factor out the Vandermonde factor in 
the denominator. 
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This method of selecting pole terms is simple but at increasing orders 
there are many terms and many coefficients. A graphical characterization of 
those terms is thus necessary, in order to extract the various contributions 
to the pole terms of Pipn- For the general problem of order k one has 

^ = e^ ''"' A(^^(Ay '^'°1 + ^™ ^^'' + ns + ■ • ■) + ^A] (^2^13 + • ■ ■) 
+ qi,i](ri2r34 + ---) + 0(r3)] (38) 

We represent the terms which are products of Tij by a graph in which one 
draws a line between the points (i) and (j) whenever a factor r^j is present. 
For a given graph, a graph-dependent coefficient characterizes its weight in 
the function /^. For instance the graph [graph]= [1,1], has a weight C[jj], 
which is the coefficient of the sum Sr„j„2''n3n45 where the Uj are all different. 
The graph [A] consists of two lines meeting at one point, and it comes with 
a coefficient C[a] in front of the sum Sr„^„2r„j„3 in which the three rii are 
different. 

One may now apply the rules (i)-(iii). 

at first order: We consider the pole term of the form . From i/j\q-] 

the operation (i) gives — 2C[o] = —2. The rule (ii) gives the same pole with 
the coefficient — 4C[i]. The rule (iii), from the terms — 2(Ai — \m)/{xi — X2) 

, where m is not equal to 1 or 2 we obtain a term —2{k — 2)C\i] . 

xi — a;2 
Therefore the cancellation yields the relation 

-2^0] - 4^1] -2(A;- 2)^1] = (39) 

which gives Cp] = — ;|- 

at order two: There are two coefficients corresponding to the two graphs of 

that order, namely C[a] and Cp^]. 

To determine CfAi, we compare the pole term AiTis or equivalently 

Xi — a;2 
1 
— AiAs^i term in ipi and ■i/j2- Applying (i) for ipi yields — 2(Ai — A2)ti3. 

Xi — X2 

Applying (ii) for ipiA], we have — 4(Ai — A2)'ri3. Applying (iii) for ■?/'[a] term, 
we get — 2(Ai — Am)ri3, where m is not equal to 1,2 or 3. Thus wefind the 
relation 

-2C[i] - 4qA] - 2(A; - 3)qA] = (40) 

which leads to C[a] = yk-i) - 

To determine Cni], we collect the pole terms of the form —2 r34. 

xi — a;2 
In T34 = (A3 — A4)(a;3 — X4), we look at the term A3X4, and thus select the 

pole term AiA3X4/(a;i — 2:2). From (i), one obtains a coefficient 2C[i] for this 

10 



pole (fig. 2 (i)). Since we are looking at terms containing AiA3a;4, they arise 
from ■?/'[A] as — 2(Ai — A4)r34 with the coefficient C[a] (fig.2(iii)). Finally from 
■^[1,1] J we have three type of contributions; the first one is Fig.2(ii),it gives 
— 4(Ai — A2)r34C[ii]. the second one is 2(A2 — A3)ri4C[ij] (fig. 2 (iv)), and the 
last one is — 2(Ai — \m)T34C[i,i] (Fig.2 (v)). In this term m should be larger 
than 4, and a factor {k — 4) arises when one sums over all possible values of 
m. 

From these terms (fig.2 (i)-(v)), we obtain the pole terms AiA3X4/(xi— X2). 
Thus adding these terms, one finds the relation 

2qi] + (4 + 2 + 2(A; - 4))qi,i] + 2C[a] = (41) 



2...1 




3 3 

(i) (ii) 

Fig.l, three contributions of pole terms from the graph C[a]. 
(i) -2(Ai - A2)ri3, (ii) -4(Ai - A2)ri3, (iii) -2(Ai - A„)ri3. 



From this equation we find 

k{k-l) 



Cm = Ju L (42) 



2 ... 1 2 1 2 ... 1 2. ... 1 2 ... 



3 4 3 4 3 '4 3/ \i 4 3 ^ 



(i) (ii) (iii) (iv) 



V 



Fig.2, five contributions of pole terms from the graph C[a] and Cpjj. 
(i) -2(Ai - A2)r34qi,i], (ii) -4(Ai - A4)r34qA], (iii) -2(Ai - A4)r34qA], 
(iv) -2(A2 - A3)ri4C[i,i], (v) -2(Ai - A„)r34. 

at third order: 

At order three , there are five different graphs contributing to ip^. 

^3 = ^[Y] + V'[N] + ^[A] + ^[A,l] + ^[1,1,1] (43) 
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The method for obtaining the coefficients C[y], C[n], C[A], C[aj], C[ijj] is the 
same as for the second order calculation. For instance, to obtain C[y], we 
look at the pole AiA3A4a;^/(a;i — X2). This singularity is obtained from the 
three graphs of fig. 3. Adding these three terms, one finds 



which leads to 



-2qA]-4qY]-2(A;-4)qY] = 



Ci 



[Y] 



k-2 



Q 



[A] 



1 



A;(A;-l)(A;-2) 



(44) 
(45) 



2...1 






_ (i) (ii) (iii) 

Fig. 3, graphs for the determination of C[y]. 



For C[7v], we concentrate over the combination AiA3a;ia;4, contained in the 
term (Ai — A2)T3ir34 for instance. The relevant graphs are 



V 7 

?>LA 3^4 




_ (i) (ii) (iii) (iv) 

Fig. 4, graphs for the determination of C[n]- 



These four graphs give -2(Ai - A2)r3ir34, -4(Ai - A2)ri3r34, -2(Ai - 
Am)Ti3r34(m > 4), and 2(A2 — A3)r3iri4, respectively. We extract the pole 
term AiA3a;ia;4/(a;i — X2) from these contributions, and the cancellation yields 



-2C[A] - 4qN] - 2(A; - 4)qN] - 2qN] = 



or 



Q 



[N] 



-Q 



1 



[A] 



(46) 



(47) 



k-i^''^' k{k-iy 

For C[A], we consider the cancellation of the pole terms with the residue 
AiA3a;ia;2. The cancellation gives 

2qA] + 4CiA] + 2{k - 3)qN] = (48) 
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which determines C[a] 

For C[Aj], the second order term gives — 2(Ai — A2)'r34r35 as residue. There- 
fore, we look at AiA3a;4a;5 for residue. There are six different graphs which 
contribute to this residue. 

Adding these terms, the cancellation gives 

-2C[A] - 4C - 2{k - 5)qA,i] - 4qA,i] + 4qN] - 4C[n] - 2C[Y] = (50) 
which determines C[a,i] 

^[^'^1 = -kik-mk-2) ^^^^ 

The last coefficient Cpjj] needs 7 graphs when we consider the cancella- 
tion of the residue AiA3A4a;5a;6. The cancellation reads 

-4C[i,i] - A{k - 2)qi,i,i] - qA,i] = (52) 

which determines C[i i^j 

C[i,i,i] = -^—j ^(-C[i,i] -2qA,i]) 

k^-Qk + 10 

(53) 



k{k-iy{k-2y 

Note that this coefficient Cpj^j has a nontrivial numerator, which is not a 
product of simple factors. This term is characterized as a non-intersecting, 
three lines graph. When the graphs are made of non- intersecting lines, the 
expression for C becomes more complex, and there is no obvious simple 
structure. In other words, for graphs with non-intersecting lines, the de- 
composition rule found in ref. [5] for n=4 does not hold. For graphs with 
connecting lines, we find that the decomposition rule does hold. 
at fourth order: 

we have computed the coefficients of every graph up to order four. How- 
ever wed have not been able to find an a priori rule giving the weight of an 
arbitrary graph. In Table A, the coefficients C are listed. 
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Table A: WKB-expansion coefRcients for /5 = 4 



1=1 


^11 = 


1 
k 




1=2 


C[A] = 


1 (^ k-2 




fc(fc-l)' ^[1-1] A:(fc-1)2 




1=3 


C[A] = 




1 


^W fc(fe-l)2' '-'[^] 


k{k-\)(k-2) 




C[A,i] = 


fc-3 /-( 


fc2-6fc+10 




- fc(fc-l)2(fc-2)' "-^W.I] - 


fc(fc-l)2(fc-2)2 


1=4 


C[n] = 


r^ r^ 1 




•-[a "-[^Al fe(fc-l)2(fe- 


-2) 




C[M] = 


A:-3 fi 


1 




fc(fc-l)2(fc-2)2' ^W fc(fe-l){A:-2)(fc-3) 




C[A,I] = 


_ r^ _ (fc-3)2 


(~i fc-4 




^[J^.i] fe(fc-l)2{fc-2)3' 


^{'^A ~ k{k-l)2{k-2){k-3) 




C[A,A] '- 


fc-4 






~ fc(fc-l)2(fc-2)2 






^[A,I,I] 


fe3-10A;2+34fe-38 
fc(fc-l)2(fc-2)3(fc-3) 






C[i,i,i,i] 


A;4-14fc3+76fc2-188fc+174 
~ fc(fe-l)2(fc-2)3(fe-3)2 





5 The r-expansion for (3= 2m, (m=2,3,4,- • •) 

Hereabove we have dealt with j3 = 4. For f3 = 2m, m=2,3,4,..., again we 
found that the WKB expansion for / terminates after a finite number of 
terms. However a given Tj^ may now appear non linearly, up to degree (m — 
1). Therefore, the graphical representations may have now multiple lines 
connecting the two points (i) and (j) (with a mulitiplicity at most equal to 
{m — 1)). Therefore there are now new types of graphs which did not occur 
for (3 = 4(m = 2). 
From the equation 

p-T.g,-P(l-i)T.T^^.-T.^l (54) 

Pip = (55) 

We expand ip as 

^ = [A(:^^A(A)]/3/2-i [^+^m(^^^+- ■ ■)+C[||](^i'2+4+- ■ ■)+qA](ri2ri3+- ■ ■)+■ • 

(56) 
Pip = connects the terms in ip of degree (n-1) to those of degree n. 
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From ipo, we have a pole term of 
From Pip [I], we obtain 

Thus we obtain the first order result, 

qi] = -i (59) 

It is remarkable that this result is the same as for [3 = 4 but now it holds 
for arbitrary p. Every residue had the same factor (| — 1) and thus P drops 
from the equation. This property holds to all orders, and this characteristics 
will become important when we discuss the duality relation for zonal poly- 
nomials. Therefore, the cancellation of the pole terms gives the same results 
as for /3 = 4, except that there are now new terms coming with the multiple 
lines. 

We list here the results of the cancellation conditions which determine 
the coefficients characterizing the r-expansion. They are valid for arbitrary 
/3, and in fact they do not depend on p. 
Second order. 

qi] + (A;-i)qi,i] + qA] = o (60) 

qi] + {k- 1)C[A] + 2qn] = (61) 

(60) is same as (41), and (61) is different from (40) by the double line term 

C[ii]- 

Third order. 

qi,i] + {k- 2)qi,i,i] + 2qA,i] = o (62) 

qii]+3qni] + (A;-l)q^] = (63) 

qA] + {k- i)C[A,i] + qY] = (64) 

qA] + (A;-2)qY]+4q^] = (65) 

qA] + ik- 2)qN] + qA] + qii,i] + q z ] = o (66) 

qA] - qii] + 2qA] - q / j + (fc - 3)qN] -{k- 3)qn,i] = o (67) 
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Fourth order. 

qi,i,i] + {k- 3)qi,i,i,i] + 3qA,i,i] = o (68) 

qA,i] + {k- 2)C[A,i,i] + qY,i] + qA,A] = o (69) 

qA,i] + {k- 3)qA,A] + 2qM] + 2qA,ii] = o (70) 



2C[ij,i] + 8C[N,i] — C[A,i] + 4C[ii,ij] 

-qA,A] - 3qY,i] + {k- 6)qA,i,i] = o (7i) 

These equations are obtained from the cancellation conditions of the 
residues of the pole ^ ^^ proportional to AiA3A4A5X6a;7a;8, A1A3A4A5XQX7, 
A^A3X4a;5X6, A^A2A3X4X5a;6, respectively. 

Further, from the residue of the form A^A3X4X5, one finds 

C[ii,i] + C'[N] + 3C[g] + 2C[nji] + C[n] 

+ {k - 2)qA,ii] + ik- 2)C[M] = (72) 

From the vanishing of the residue containing A^A3A4a:5, one obtains 

C[Y] + 2q>] - 2Ciui + 2q /^i] + 2q^] + ik- 3)^//] = (73) 

and we have 

C[Z,i]=^[ui (74) 

However, starting with order three and higher, the coefficients of the r 
expansion are not determined uniquely, since the r variables are not inde- 
pendent. Starting with degree three there are identities, i.e. polynomials 
vanishing identically. In the appendix C, we discuss the origin of these iden- 
tities and derive the cubic identity and the quartic identities. 
The cubic one is 

/3 = [II,I]-[N] + (k-3)[A]=0. (75) 

and the quartic ones are 

2[II,I,I]-[N,I] + (k-5)[A,I] = (76) 

and 
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[III,I] +2[II,II] - [>] + (k-3)[ A] -4[n] - [U] = (77) 

and 



[A, II] - [M] + 2[A,I] + (k-4)[n] - l(k-4)[II,II] = (78) 



and 



[Z,I]-2[ZZ]-2[A,I] + (k-4)[>]-2(k-4)[n] + (k-4)[II,II] = (79) 

Those identities give some freedom in writing the above expansion. In- 
deed one may add one of those identities multiplied by an arbitrary coefficient 
to the r expansion of the HIZ-integral. For the third order, the coefficients 
Cii,i,C'n and Ca may be shifted as 



ClI,I — ^ C*!!,! + OL 

Cn — > Cn — ol 

Ca -^ Ca + (A;-3)« (80) 

where a is arbitrary. 

At order four, from (76), (77), (78) and (79), the following shifts leave the 
expansion unchanged 



Cii,i,i - 


-^ Cn,i,i + 274 


Cn,i - 


^ Cn,i - 74 


C*!!!,! - 


^ Ciii,i + 7i 


Cii.ii - 


-^ Cii,ii + 271 - - (A; - 4)72 + (/c - 4)73 


c> - 


-^ C> - 71 + (/c - 4)73 


Ca - 


-. Ca + (A;-3)7i 


Cn - 


-^ Cn - 471 + (A; - 4)72 - 2{k - 4)73 


Ca,i - 


- Ca,i + 2(A; - 4)72 -2(A;- 4)73 + (A; -5)74 


Cu - 


-^ Cn - 7i 


Cm - 


-* Cm - 72 


cu - 


-^ C/^i + 73 


Cll - 


-^ C//-73 



where 71,72 and 73 are arbitrary numbers, arising as multiplicative factors 
of the quartic identities. 
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In Table B, we give the coefficients C[graphs]- For the coefficients, which 
have ambiguities due to the above identities, we use a fixing rule which will 
be discussed in the next section. 

The values of Table B are also derived systematically by the r-trasformation 
of the Jack polynomial series as discussed later in Appendix B and Appendix 
D. In Table B, we use the parameter a instead of f3. The parameter a is a 
parameter of Jack polynomial, and in Table B, a is given by a = ^^. For 
P = 4, a becomes —1, and the values in Table B coincide with the values in 
Table A. 



Table B: WKB-expansion coefRcients for P = 2m 



1=1 


qi] = 


1 

k 


1=2 


C[A] = 


1 


k(k+a) 




C[i,i] = 


^ (I \ "1 
" k{k+a)^^ ' k-l> 




C[ii] = 


1+a 




2k(k+a) 


1=3 


C[Y] = 


1 


k(k+a){k+2a) 




C[A,I] = 


1 fl 1 2a ^ 
k{k+a)(k+2a)y'- ' k-l) 




<^[III] = 


(l+a)(l+2a) 
f!,k{k+a){k+2a) 




'^[l]~- 


l+a 




2k{k+a){k+2a) 




C[i,i,i] '- 


1 /I 1 3a 1 2a2 x 




k{k+a){k+2a)\^ ' fc-1 ' (fc-l)(fc-2)/ 




Qll,l] = 


1+a /-, 1 2a N 
2fc(fc+a)(fc+2a) V-^ ' k-l> 




C[N] = 


A:(A:+a)(fc+2a)V-^ ' k-l' 




C[A] = 


1 fl "M 

k(k+a){k+2a)\^ k-l) 



(The last three coefficients C have arbitrariness due to the cubic identity) 
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1 A/^ 1 

i— "±<-^[X] — k(k+a)(k+2a)(k+3a) 

^[Y,I] = k{k+a){k+2a){k+3a) (^ + k^) 

^[A,A] = k{k+a){k+2a){k+3a)y'^ + fe+a-J^^ + fc^) 

/^ 1 /I 2a I 7a i 6a^ 2a 

•-"[A,!,!] — fc('fe+a)(fc+2a)(fc+3a) V-"- fc+a-1 "•" fe-1 "•" 



fc(fe+a)(fc+2a)(fc+3a) V fc+a-1 ' fe-1 ' (fc-l)(fc-2) (fc-2)(fc+a-l) - 

/^ 1 /I _|_ 6a J 9af i 8af 

'-'[1,1,1,1] ~ A:(A:+o)(fc+2a)(fc+3a) V "*" k+a-1 "*" (fc-l){A:+a-l) "•" (fc-2)(fc+a-l) 
I 25a3 8a3 , Ga" \ 



(fe-l){A:-3)(fc+a-l) (fc+a-l)(fc-2)(A;-3) ^ (fc-l)(fc-2)(fc-3)(A:+a-l) ^ 

1+a I'l , 2a 



/^ 1+a /'I _|_ _2c^\ 

^E] ~ 2A:(fc+a)(fc+2a)(fc+3a) V "•" k-l^ 

r = (l+a)2 

^[<] 4fc(fc+a)(A:+2a)(A:+3a) 

/^ 1+a 

^[H ~ 2fc(fe+a){fe+2a){A:+3a) 

^ ^ (l+a)(l+2a) 

^ [ i. ] 6fc(fe+a){A:+2a)(A:+3a) 

^ ^ (l+a)(l+2a)(l+3a) 

^[IIII] 24A:(A:+a)(fc+2a)(fc+3a) 

(These 10 coefficients are determined without ambiguities) 
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^[Hl ~ 2fc(fc+a)(fe+2a)(fc+3a) V "•" k-l ' 



_ l+o /I _ 2af_^ 

[^] ~ 2fc(fc+a)(fc+2a)(fc+3a) V k-l' 



^\LL\ fc(fc+Q)(fc+2a)(fc+3a) ^ "*" fe-l/ 

/^^ , _ 1 {^ _ "("^1) '| 

'-'ly fc(fc+a)(A;+2a)(A;+3a) V-^ fc-1 / 

[i..I] ~ 2fe(A:+a)(fc+2a)(fc+3a) v-^ "•" k-\l 

(~i l+g /I I a 'i /-r i 3a \ 

'-"[A.II] — 2fc(A:+a)(fc+2a)(fc+3a) V-^ """ fc+a-l/V-^ "•" k-\l 



1+a /I 2a I 7q I 6a^ 2a^ 



/^ ^ l+g /I 2a _|_ la i 

'-'[11:1.1] ~ 2fe(fc+a)(fc+2a)(fc+3a) V-^ fe+a-l """ fc-1 """ (fe-l)(fc-2) (fc-2)(fc+a-l) ' 

fi ^ {\+af /^ I J^ _| 2a^(2+a) x 

'-'[II.II] 4fe(fc+a)(fc+2a)(fc+3a) V "*" fc-1 "*" (l+a)(fc-l)(fc+a-l) '' 



fi ^ (l+a)(l+2a) /-, , 3a \ 

^[III.I] 6fc(fc+a)(fc+2a)(fc+3a) V-"- ^ fc-1/ 

(~i __ 1 /I _|_ 2a \ 

^[°1 ~ fc(fc+a)(fc+2a)(A;+3a) V """ fc-1'' 

^[M] = fc(fc+a)(fc+2a)(*;+3a)''-'- + fc+a-Jv-^ + fcTT) 

^ _ 1 /I I _8^ _ Jo of I a(l+a)(4fc+3a-4) ■ 

^[N,I] fe(fc+a)(fc+2a)(fc+3a) V-"- "•" fc-1 fc-2 (fc+a-l)(fc-2) "•" (fc-l)(fc-2){fc+a-l) • 

^[^'I] ^ fcffc+a)ffc+2a)ffc+3a) ^ 



1 /-, a(a-3) a^(4a-3) a^(3a^-2a+3) 



fc(fc+a)(fc+2a)(fc+3a) V fc+a-1 (fc-2)(fc+a-l) (fc-l)(fe-2)(fc+a-l) ' 

(These 13 coefficients involve arbitrariness due to 4 quartic identities. 
The determinations of these coefficients are given in Appendix D.) 



\L\ [E] [A] [H [Z,I] 

\LL\ [>] 

Fig. 5. symbols for the several graphs. 
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6 The r-expansion from zonal polynomials for 

The Table B shows the first coefficients of the WKB-expansion for general 13. 
It was derived pertubatively assuming that (3 was an even integer. However 
the expression may be analytically continued to all integers. In this section 
we consider f3 = 1, which is of practical importance since it is the case of a 
measure invariant under the orthogonal group (the matrices X and A being 
real and symmetric). 

The HIZ integral (1) may be expanded in a series involving products of 
zonal polynomials Zp{X) and Zp{A), as 



'/3=1 



[exptr{XgAg ^)dg 

OO 1 „ 

^ — / [tr(XgAg-i)]-dg 
_n m\ Jo(k) 



m=0 ""• •'Oil^'l 



E A^ E Xp(l) / ZpiXgAg-')dg 
f— ^ m! 2m ! , , Jok) 



m=o^'(2m)!;^) J Oik) 

where p{m) is the partition of order m ; Zp{X), the zonal polynomial, is a 
symmetric homogeneous polynomial of degree p in the k eigenvalues of X 
[15]. The third equality in (82) follows from the identity 

{trMy = c,J2Xpil)Zp{M) (83) 

p 

where Cg = J2 Xp(1)- The sum over p runs over the partitions of the number 
q, i.e. over all Young tableaux with q boxes. For the orthogonal group Cq = 
2'^q\/{2q)\. The function Zp{M) is a symmetric function of the eigenvalues of 
M of degree q. 

In the unitary case, a similar expression in terms of a character expansion 
is well known ; it is used explicitly in [12] ; there Xp(1) is the dimension of 
the representation of the permutation group of p objects corresponding to a 
given Young tableau. (In an appendix, we give the general expression of this 
integral for arbitary j3, including the unitary case.) 

The zonal polynomial Zp{X) beeing an homogeneous symmetric function 
of xi, ...,Xk, may be expressed in terms of the sums s„ (n=l,2,3,...), 

s„ = trX" = x'l + x^ + ■ ■ ■ + x"^. (84) 
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The zonal polynomial Zp{X), the constants Xp(l)) ^-iid Zp{l), which is ob- 
tained by setting all Xi = 1, are hsted up to order five in the Appendix A, 
when one takes a = I = 2. 

From (21) and (82), we have the relation, 

eS-i"'^7/3=i = //3=i(a:„A,) (85) 

For /9=l,the Vandermonde factor disappears, since it is raised to the power 
is j3 — 1. For arbitrary (3 there is a sum in (1) over the k\ permutations of the 
Aj. However, if one expands in powers both the exponential term, together 
with //3=i, it turns out that each term of given order is a symmetric function 
of the Xi and of the Xi. Therefore the sum over permutations in (21) is not 
necessary for (3 = 1. 

In the previous section, using the coefficients C of the Table B, setting 
f3 = l,we have obtained the r-expansion for ffs=i. The same result may be 
derived from the zonal polynomial expansion for (3=1. 



Indeed, using the explicit values of the character Xp(l),the zonal polyno- 
mial expansion (82) reads 

\2 I o„ i^wii^ /'\^^2 



h=i 



+ 2 



s,{x)s,{\) l Msr{x)f + 2s,{xms,{\)f + 2s,{\)) 
k 6^ k{k + 2) 

{{s,{x)f - s,{xms,{\)f - s,{\))^ 



k{k-l) 



^ 1 r Z^3](x)Z[3](A) ^ ^ Zi2i]{x)Z[2i]i\) ^ ^Z^i3](a:)Z[i3](A) 



90'A;(A; + 2)(A; + 4) A;(A; + 2)(A; - 1) A;(A; - 1)(A; - 2)' 
+ ■■■ (87) 

For a comparison of this expansion with our earlier expressions ff3=i , one 
must still expand the exponential factor e^^*^'. But this exponential factor 
has no explicit k dependence ; thus the k -dependent coefficients in /^=i 
will not be modified when one expands the exponential. For instance, at 
first order, the coefficient | is present both in (86) and (87) ; expanding the 
exponential e^^^^% we do find that the two expansions coincide. At order 
two, one writes the coefficient of T12T34 as 

fc + 1 - ^ ^ I ^ ^ r88^ 

A;(A; + 2)(A;-1) 3k^k + 2 k-l'' ^ ' 
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then the terms of degree two in e^^^'^^f/3=i have coefficients which are either 
kd+n *-'^ k(k-i) ' ^^^y ^^^ exacly the inverse of the dimensional constant 
Zip]{l) in (87^. 

At order three, consider for instance Ti2T34T^q (three non-intersecting 
hnes). Taking in Table B the coefficient of this term, one may decompose it 
as 

k'^ + 3k-2 1 9 

+ 



A;(A; + 2)(A; + 4)(A;-l)(A;-2) 15A;(A; + 2)(A; + 4) 15A;(A; + 2)(A; - 1) 

+ - f89) 

15A;(A;-l)(A;-2) ^ ' 



which is the sum of 



9 5 



15^[3](4) 15^[2i](4) 15^3] (4)' 

The numerators of these terms are in the ratios 1:9:5, which are also the 
ratios of the characters X[3](l) • X[2i](l) : X[i3](l)- (We will soon find the 
reason for this coincidence.) 

Let us return to the derivation of the r-expansion from zonal polynomials. 
We ffist shift the diagonal matrices X and A to make them traceless, 

(90) 
(91) 

(92) 
Noting that tr A = tr X = 0, the HIZ integral is 









A 


= -trA + A 
k 








X 


= ^trX + X 
k 


In terms of 


eigenvalues. 


it 


is 


lU^a-Xb 



From the expressions of zonal polynomials in terms of the symmetric 
functions s„ given in Appendix A, noting that Si = J2xa = 0, and using 
simple identities such as ^(E Xi){T,Xi) = E AjXj - ^ Ei<j(Ai - Xj){xi - Xj), 
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the integral becomes 



/ = eS^-iE.<.-[i + ^__^,,(,),,(A) 



{k + 2){k 

2k{e + k + 2) ~ 

^ {k + l){k + 2){k + A){k + Q){k-l){k-2){k-?,f^^'''^^^ ' 

k^ + 5k^-6P-36k + 72 -^ 

^ 2k{k + l){k + 2){k + 4){k + 6){k -l){k- 2){k - 3) ^^2 WS2(A)) 
2(2A;2 + 3A; - 6) 
(A; + 1)(A; + 2)(A; + 4){k + 6){k - l){k - 2){k - 3) 
X {{s2{S:)fs,{\) + s,{S:){s2{\)f) + ix% (94) 

The paired product S2{x)s2{X) is expressed in terms of the Tij as 

rCt ■ , ■ rCt ■ , ■ 

(A; -1)2 2 2(A;-1)^ 2^ 

= ,2 -L^^j p E^u^ifc + pE^^j-^H (95) 

Aj ■ , ■ Aj IXi 

The second sum is restricted to j < k and i,j, A; all different. The last 
sum is restricted to i < j,k < I and i,j,k,l all different. The identity 
(95) holds for arbitrary Xi and Aj. In order to fix the coefficients of this 
identity, one may choose simple values of X and A. Let us, for instance, take 
X = A = (1, 1, ..., 1, 0, ...0), where the eigenvalue one is q-fold degenerate, 
and zero {k — q). Then we find 

S2{X)S2{X) = ^^^^^ (96) 

Note that — is invariant under the replacement q —^ k — q. The 

k^ 
three terms of the products of r in the r.h.s. of (95) become q{k — g), 

q{k — q){k — 2)/2 and q{q — l){k — q){k — q — l)/2 respectively. These numbers 

are also invariant by g ^ A; — g. Then with three unknown coefficients, a,b 

and c, we write 

''^\,'^^'=a[qik-q)]+,f^^-'^f-^\+cf^^-'^^^-f^^-^-\ 

(97) 

This relation faxes the three unknown coeffacients : a = — , b = — 

A;2 k^ 

k^ 
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In the Appendix B, we derive this identity with the help of differential 
operators. In the same appendix, we discuss the general expression of product 
of symmetric functions s„(x)s„(A) in terms of r. 

Remarkably, as may be checked in the explicit expressions hereabove for 
the second order, although each of the three quadratic expressions in r in 
(95) is not invariant under the permutations of Aj, the combination of the 
three types of products of Tij is a totally symmetric function as it should. To 
study further this invariance under permutation, we consider 

X = (i,...,i,o, ,0) 

A = (0, ,0,1,. ..,1) (98) 

where the eigenvalue one is q-fold degenerate, and we choose q < k. This is 
a particular permutation of the previous choice of X = A = (1, ...1, 0, ..., 0). 
We assume q < k/2. The various r are not the same as in (97) ; the three 
quadratic terms are now respectively [g^], [g^(g — 1)] and [g^(g — 1)^/2]. For 
this particular permutation of the Aj, the identity 

is indeed consistent with (95). 

The equation (99) is interesting since, in this choice of permutation, the 
value of the r's are all function of q and not of k.If we set k=q, this equation 
shows a non-trivial cancellation for three terms, since the left hand side 
vanishes ; in this case the right hand side of ((99) vanishes as 

(l-2 + l)g2(g_i)2 = (100) 

To summarize we have found a method to write S2(5)s2(A) in terms of r's. 
Given that it is a function of the r's the method consists of (i) evaluation of 
S2{x)s2{\) for the choice X = (1, ..., 1, 0, ..., 0), A = (0, ..., 0, 1, ..., 1), where 1 
is q-fold degenerate. (ii) Decompose this value as I2 [(A; — 1)^ — 2(g — 1)(A; — 1)-|- 
(g — 1)^] . (iii) Evaluate the r terms for this same choice as a function of q. (iv) 
The comparison with the expansion of (iii) fixes the unknown coefficients. 

Similar identities expressing the symmetric functions s„(X)sn(A) in terms 
of r's hold at higher order. At order three since Si{X) = 0, Zp{X) is given 
only by s^lX) (see Appendix A). In order to express 53(0;) S3 (A), into the 
Tij = {xi — Xj){Xi — Xj) polynomials, in the Appendix B, we present a method 
for deriving these identities based on a differential operator D\'l^.^. Let us 
for instance quote the result for the third order 

s,{X)s,{K) = -^^{k-lf{k-2)\rl, + ...) 
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(k-l){k-2)\ , 

(k — 2)^ 
+ 3^-^(ri2r23ri3 + ---) 

^^ {k-l){k-2) ^ 

- 12 — (ri2ri3ri4 

(k — 2)^ 

- 6 — — — (ri2r23r34H ) 

+ 12^-^(ri2r23r45 + ---) 
24 

- Tj(T12T34T56 + ■ ■ ■)■ 



(101) 



In this expression we have used the freedom given by the cubic identity among 
the r's to remove the ambiguities. 

However the result may be obtained also by the direct method which was 
used above for the quadratic case. One first choose X = A = (1, .., 1, 0, ..., 0) 
with eignevlue one q-fold degenerate ; then S3(X)s3(A) = g^(g — A;)^(2g — 
k)^/k^, which is invariant by the substitution oi q -^ k — q. For instance, 
the last term of non-intersecting graph of(ri2r34r56 + perm.) becomes q{q — 
l){q — 2){k — q){k — q ~ l){k — q — 2)/6. Each of the 8 terms is invariant by 
the substitution oi q ^ k — q. Therefore the coefficients have to be functions 
of k only. Then, we apply again the permutation A = (0, ..., 0, 1, ..., 1). The 
various r terms are polynomial in q. The remarkable factorization identity. 



k^ k^'Q 

^4 [g^(g-i)%-2)]-L ^, 



12^[g^(g - mq - 2)] - 9^-^-^[q\q - if] 



- ''-''J-''\ .^^ (102) 

extension of the analogous quartic identity, fixes the unknown coefficients. 
As discussed in the previous section, the cubic identity 

^3 = [^iV34] - [^2^-23 7-34] + {k- 3)[ri2ri3r24] = (103) 

where [t^2''"34] ~ '''f2''"34 + ' ' ' leads to an ambiguity, since one may add 13 with 
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an arbitrary coefficient a. At order three, one has from (94), 

f (3) ^ \ 

6A;(A; + 2)(A; + 4)(A;-l)(A;-2) 

X \{k^ + 3A; - 2)g^ + (-12A; - 24)g^ + (eA;^ + Vlk + 56)g^ - 48g^ + %k^q^\ 

(104) 

By taking the choice X = (1, ..., 1, 0, ...., 0) and A = (0, ..., 0, 1, ..., 1), we 
evaluate [Ti2r34r56] = |g^(g — l)^{q — 2)^. By comparing the highest order of 
q, q^, we find 

r^ - k' + 3k-2 

^P'^'^l " ~kik + 2)ik + A)ik-l){k-2) ^^^^' 

Subtracting this Cpj jj [ri2r34r56] from l'^^\ we have 

?{3) ^ I 

«' 6A;(A; + 2)(A; + 4)(A;-l)(A;-2) 

X [{6k + 18)q^ + {-7k - 41)g^ + (12A; + 12)q^ + {Ak - 4:)q^]{106) 

The quantity [ri2ri3r45] becomes g^(g — iy{q — 2) and it is order g^. Thus 
we determine its coefficient C[ai] from the q^ term in (106), 

^[^'^1 - ~k{k+m+m^ ^ ^ 

We subtract C[aj] [ri2ri3r45] from (106), and we obtain 

f(3) ^ {I7k + 31)g^ + (-18A: - 78)?^ + (16A: + 32)q^ 
"' 6k{k + 2){k + 4){k-l) ^ ' 

The terms, which give order g'', are [t^2''"34] = ?^(?^1)^) [T12T13T34] = q^{q—^Y 
and [Ti2ri3ri4] = \q^{q — l)(g — 2). ([ri2r23ri3] becomes vanishing.) Thus we 
have 

-^Y] + qii,i] + qN] = -6^(^ + 2)(A: + 4)(A:-l) ^^°^^ 

Since C[y] is 

^[^5 = "A;(A: + 2)(A; + 4) ^^^^^ 

we have 

^^"•^1 + ^t^l = "2A:(A: + 2)(A^+4)(A:-1) ^^^^^ 



We write 

5A; + 11 9 



2A;(A; + 2)(A; + 4)(A;-1) 10A;(A; + 2)(A; + 4) 5A;(A; + 2)(A; - 1) 
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In principle, the coefficients C[u,i], C[n] and C[a] have ambiguities, since there 
is an identity of (103). One method of fixing this ambiguity may be large k 
assumption: We assume 



^™ " ~2^k{k + 2){k + 4) + A;(A; + 2)(A;-1)] ^^^^^ 



and 



^[^]""[a;(A; + 2)(A; + 4) + A;(A; + 2)(A;-1)] ^^^^^ 

We assume that C[u,i] ~ |p- in the large k limit The factor | is the degereracy 
factor = (/?/4 — l)/{j3/2 — 1) for /? = 1 ( in the next section, we will discuss 
this factor by the substitution of a = 2/(2 — P)) . And C[n] ~ -p- also. Then, 
the comparison with (112) gives 

3 9 

-a + c= — (115) 

2 10 ^ ' 

Still one parameter a or c remains undetermined. The coefficient C[a] is de- 
termined by the residual equation of (113). However, the ambiguity constant 
a remains. When we consider the case f3 = 4, we have a definite value of 
C[A], since there is no double line in /9 = 4, and no cubic equation exists. For 
the general value of C[a], extending the case of /9 = 4, we write 

1 a^ 

^[^1 = ~kik + a){k + 2a)^^ ' I^] ^^^^^ 

Then, using the residual equation of (66), and (67), we determine the con- 
stant a and c in (113) and (114) as 

a = i, c=^ (117) 

which also satisfy the large k limit condition (115). In other words, we 
have assumed the value of C[a] as (116), since we are free to choose of the 
coefficients of the cubic identity. 

Note that when we subtract [C[n,i] + C'[n]]?^(? — 1)^ from (108), we obtain 

A, (3) 2g^ + 12g3 + g^ 

where there is no -j^ factor. This quantity coincides with the sum of 
C[iii][^i2]) C*!/] [^12^13] and C[y] [7-12X13X14], since they are 

5 1 

^™ = "2A:(A: + 2)(A: + 4) ^^^^^ 
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^[^' = "A:(A: + 2)(A: + 4) ^^^^^ 

where [rf2] = g^[rf2^34] = 2g=^(g - 1) and \ti2Ti^tu] = \q^{q - l)(g - 2). 

It may be interesting to investigate the fourth order term. In (94), we use 
the choice that X = A = (1, ..., 1, 0, ..., 0) where 1 enries q-times, entries 

(k-q) times. Then, noting that S2{x)s2{X) = t^ , 

,~, r., q\q-kn2q-kf ~ g^ ^k - q^k' - 3kq + Sq^ 
S3{x)s3{X) = , S4(a:)s4(A) = ^ , 

[S2[X)S2[\)\ = -^ , 

[S2(£))^S4(A)+S4(£)(S2(A)) 

_1 V-> q 

we have the fourth term together with the expansion of e *= ^ 



~NA2 /xA f~^f /xxx2_ V(^-g) (^ -3A;g + 3g2 



A;5 



/(4) 



1 



24A;(A; + 1)(A; + 2){k + 4)(A; + Q){k - l){k - 2){k - 3) 
X [{k^ + 7P + k^ - 35k - 6)g^ + (-24P - 144A;^ - 72k + 144)g^ 
+ (12A;^ + 72P + 308A;2 + 1064A; + 864)g6 
+ (-240A;^ - 1008A;2 - 2304A; - 1728)g^ 
+ (44A;^ + 188A;^ + 1944A;2 + 2064A; + 96)g^ 
+ {-576k^ - 672k^ - 192k)q^ + (48A;^ + 48P + 96k^)q^] (122) 

There are 23 terms in the fourth order for the r expansion as shown 
in Table B. For the choice of X = (1, ..., 1, 0, ..., 0), A = (0, ..., 0, 1, ..., 1), 
each r terms are evaluated as they are polynomials of q, and each terms 
have different orders of q. The highest order of q is obtained from the term 
[T12T34T56T78] , (here [r ■ ■ ■] means the sum of the permutation of the indecies). 
It becomes 

[mmnem] = ^g'(g - l)'(g - 2fiq - 3f (123) 

The order (f term is only this term. Then, from (122), we obtain the coeffi- 
cient of this term \r\2T3i^T^(,T'j^ , as 

^ (A: + 3)(A:^ + 6A;+1) 

''''''''' ~ k{k + l)(fc + 2){k + 4)(A; + 6)(A; -\){k- 3) 

1 f ^2 ^"^ ^ ^ 



A;(A; + 2)(A; + 4)(A; + 6)^ k-\ {k-\){k-'i) (A; + 1)(A; - 3)^ 

(124) 
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This result is expressed by the use of the dimensional constant Zp{T) and by 
the constant Xp(l)) which appear in the zonal polynomial expansion, 

The next highest order term of g^ is [ti2Ti3T45Tqj], which has a factor 
C[Axi]- The value of [ti2TisT45Tqj] is |g^(g — l)^(g — 2)^(g — 3). Then from 
the coefficient of g^ in (122), we obtain the coefficient C[a,i,i], 

-0[A,i,i] - -0[i,i,i,i] 



2-,M,ii 2^^'^'^'^J 2Ak{k + l){k + 2){k + A){k + Q)ik - l){k - 2){k - 3) 
X (-24A;^ - 144A;2 - 72A; + 144) (126) 



which reads 



k^ + SP + 13A; - 2 



k{k + 1)(A; + 2){k + 4){k + Q){k -l){k- 2) 
1 r 10 20 



k{k + 2){k + A){k + 6)' k-l {k-l){k-2) {k + l){k-2)' 

(127) 

We note here that the following relation of (68), which was derived for 
the arbitrary value of (3 in the section 5, is satisfied by above two expressions, 

qi,i,i] + 3qA,i,i] + ik- 3)qi,i,i,i] = o (128) 

since we have from Table B, {P = 1), Cpjj] = - k(k+2)(k+4) (^ + fc^ + 
The pole term ^^ in Cpjjjj is cancelled by the pole term of 



(A:-l)(fc-2)/- -^^^"^ t'^'^^ uv^xiii ^_|_j^ 

C'[A,i,i]. This -^ factor comes only from the dimension of the zonal polyno- 
mial Z[22], and it appears ffist in the fourth order terms. 

For the terms, which give order q^, new 4 terms [Ti2ri3r34r5g] = q'^{q — 
l)2(g-2)2, [n^risTMrse] = \q\q-l)\q-2){q-?,), [T^2mu^T,,] = \q\q- 
lf{q - 2f{2q - 5), [(^2)^34X56] = \q\q - lf{q - 2f appear. 

From (122), after the subtraction of the contribution of [ti2T^at^qT'j^] and 
[''"i2Ti3T45r67], the q^ part of (122) becomes 



31A;3 + 196A;2 + 119A;-310 
12A;(A; + 1)(A; + 2)(A; + 4)(A; + 6)(fc - l)(fc - 2) 



Cg6 = ^^,^^,^ ^ ^^^,_ ^ ^^^,^ ^ ^"^^,^ ^ ^^^,_ :^^ — ^ (129) 



which have to be the sum of these four terms, to determine these four coef- 
ficients, we fist fix two coefficients C[y,i] and C[a,a] by the residual equation 
of (69). 
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These are satisfied by the values of 

^[^•^1 = A;(A; + 2)(A: + 4)(A: + 6)^^ + I^^ ^^^"^ 

1 8 8 

^f^'^l = A;(A; + 2)(A; + 4)(fc + 6)^^ + I^ + (A; - 1)(A; + 1)^ ^^^^^ 

Then, we subract these two contributions from (129), and the remaining 
two coefficients are 

C[N,i] + 2^[ii,i,i] - 4A;(fc + l)(fc + 2)(A; + 4)(A; + 6)(A; - l){k - 2) ^^^^^ 

From the residual equation of (71), we are able to confirm above result. 

There are ambiguities in the fourth order by the 4 quartic identities as 
shown in appendix C ((282) ~ (285)). Due to the quartic identity (282), it 
is free to choose the value of Cpi j j] . Here we assume that C[n j j] is same as 
C[A,i,i] except a factor of the multiple line, |. 

^ ^A A:3 + 8P + 13A:-2 

^"''•'^ ^2'k{k + l){k + 2){k + A){k + Q){k-l){k-2) ^ ' 

Then we have 

A;3 + 6A;2 + 3A; - 6 
k{k + 2){k + A){k + Q){k + l){k - l){k - 2) 

These two values satisfy (132). 

The limit of g -^ 1 in (122) determines uniquely the coefficient of the 
term (ti2)'^, which is four line degeneracy, and it has a value q^. Other terms 
are proportional to g — 1. The coefficient, therefore, becomes 

35 

lk{k + 2){k + A){k + Q)' 

The limit g ^ 2, in (122), after the subtractions of [ti2]^, becomes 

f (4) ^ 125P + 320A: + 291 

2k{k + l){k + 2){k + 4){k + 6){k - 1) ^ ' 

whichisthesumof [ri2r23ri4r34], [ri2(ri3)V34], [(ri2)V23r34], [{^2^^^], [{Tufris], 
[(ri2)V34],[(ri2)^(r34)^], which have values of lq^{q - l)^g^(g - l)^2g2(g _ 
l)^,g^(g — l),2g^(g — l),g^(g — 1)^, and |g^(g — 1)^, respectively. By the sum- 
mation of these terms with the coefficients we have for q=2, 

/(') = TT, o^// .^n ^[(1 + T^) + 6(1 



Cm] = ,,, , ^,,, , ,J, , ^J, , ,,,, ,,,, ., (134) 



Cpiii] = 7^771— 1W71—^V71—1^- (135) 



A;(A; + 2)(A; + 4)(A; + 6)'' k - l' ' k-V 

4 , , 6 , 9, 8 32 

+ ^2(l + ^) + 9 + 20 + 10(l + ^) + -(l + ^ + ^^^-^^^^^)] 

(137) 
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where 7 terms are added respectively. We derive the following expression, 
which involves -j^ pole. Other terms are easily derived. 

^["'"]"4A;(A; + 2)(A; + 4)(A; + 6)^^A;-l + 3(A;-l)(A;+l)^ ^'^^^ 

and 

5 6 

^P"'^l = 2A:(A: + 2)(A: + 4)(A; + 6)^^ + ^^^ ^^^^^ 

We have following equations, 

C[iii] + 2q ^ ] - q A ] + C[ u ] + (A; - 2)qni,i] = (140) 

q Z ] - C[iii] + (A; - 3)q u ] - (A; - 3)qni,i] + 2q a ] = (141) 

One could eliminate the coefiicientof C[/\], the triangle diagram with one 
double line. Then, we know that C[in,i] has no pole of j^. This is consistent 
with the expression (139). 

We have up to now confirmed 14 terms at fourth order. The number 
of fourth order graphs is 23, and among them, three graphs exist, which 
contain the triangles. These terms are vanishing for the present choice of 
X and A. Thus, 6 terms still have not yet been determined. These 6 
terms are [ti2Ti3TuTi5],[ti2Ti3TuT45], [ri2r23r34r45],[(ri2)Vi3ri4],[(ri2)Vi3r45] 
and [7-12X13 (r45)^]. They have values of ^g^(g — l)(g — 2)(g — 3), q'^{q — l)^{q — 
2),q\q-lf{q-2),q\q-l){q-2),2q\q-lf{q-2),^nd2q\q-l)\q-2). 
The sum of these 6 terms should be the value of (122), once one subtracts 
the subtractions the previous 14 terms. 

From the equations of the residues, we obtain 

r ^r {k + 3){5k + 2i) 

^[^'"1 + ^[^1 - 2A;(A; + l)(A; + 2)(A: + 4)(A; + 6) ^^^^^ 

which satisfy the residue equation of (70). Note that C[m] is not uniquely 
determined in views of the quartic identity discussed in appendix C 

For fixing this ambiguity, we take the value of C[a,ii] to be the same as 
C[A,A] within a factor |, 



k^ + eP-k- 3t 
2 k{k + l){k + 2){k + A){k + 6) 



C[A,ii] = 7^ ,,,,, , ,,,T , ^,„, , ,,„, , ^, (143) 



Note there is no -j^ factor in the expression as same as (142). From (142), 
C[M] is determined as 

r = (fe + 3)^ n44^ 

'""^ A;(A; + 2)(A; + 4)(A; + 6)(A; + 1)(A;-1) ^ ' 
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In the next section and in Appendix D, we discuss other terms for general 
P (or a). 

As a summary of this section, from the zonal polynomial expansion of the 
HIZ-integral, we have obtained the series given by the symmetric polynomials 
of (94). By writing these symmetric polynomials s„(x) in terms of r, we have 
obtained the r-expansion of the HIZ-integral in the form e^^*^'/, where / 
is given as (86). We have found that the coefficients are expressed by sums 
of the inverse of the dimensional constants. We have shown that e^^*'*'*/ 
is invariant under the permutations of Aj. Therefore, the summation of the 
permutation contained in (8) is not necessary for (3 = 1. 

7 Character expansion for general /3 and a 
duality 

In this section, we generalize the results of the previous section to arbitrary j3, 
and we discuss the derivations of the r-expansion using a dual representation. 
We write the zonal (Jack) polynomial expansion for the integral Ij3, 

where we have shifted Xi —* Xi = Xi — ■^ Y^i=i ^i (subtraction of the mean). 

The generalized zonal polynomial Zp{X) has a parameter a, which is 
equal to the present |. This polynomial is called as Jack polynomial. As for 
/5 = 1, we expand in a power series for small X and A. The quantity, which 
is an generalization of the integral (1) to arbitrary f3, is a symmetric function 
by interchange of the Xi and Aj. We shall write this quantity /^ in terms of 

Ip = eEti->A,[^ ^ ^^^^[^^j ^ /^[I,i][ri2r34] + ■ • ■] (146) 

This expansion is different from the form discussed in section 2, which 
was 

h=Y: ^u , ^'t_i [1 + C[i] [^12] + qi,i] [ri2r34] + ■ ■ ■] (147) 

■perm. lli<jVJj/ 



The differences are that (147) has a Vandermonde determinant, and also it 
has a sum about the permutations of the Aj. The expression (146) has no 
Vandermonde term and no sum over permutations. Therefore, the coefficients 
C[3rap/i] and Digraph] are not the same in these two expansion. (The P = I 
case is exceptional since there the two expansions coincide). 

To obtain the coefficients D in the expansion (146), we take X and A as 
X = (1, ..., 1, 0, ..., 0) and A = (0, ..., 0, 1, ..., 1), where 1 is q-fold degenerate, 
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and q < k/2. Since J2i<jTij becomes —q^, with this choice, we have up to 
second order, 

/ = ei:-^ii + ^ + ^ + ...] 
^ f' + 2(rT^(^[^i(') z„(i) +^'^'^'^ — viw — ^'■■■■^ 

The symmetric functions s„(x) = I]a:" are easily evaluated within this 

choice of X, and we obtain S2(x) = g(l — — )^ + (A; — glf— ^)^ = ; , 

k k k 

1 , ,~\ /AN q^{k-qY 

and thus S2(a;)s2(A) = -7. . 

k^ 
Thus we have 

We denote Z^perm. (''"12 + ■ ■ ■) by [ti2]- Under this choice of X and A, we 
find that [rijTki ■ ■ ■] are polynomials in q. 

We evaluate them at order two as [ri2r34] = \q^{q—lY, [T12T13] = 9^(9— 1), 
[(''"12)^] = (f- From (146) and (148), we write 

1 + £>[!] [ri2] + /^[i,i] [Ti2T^i] + Z^[A] [ri2ri3] + /^[ii] [(ri2)^] H 



1 + Dyr^{-q^) + /^[i,i](ig'(g - 1)') + /^[A]g'(g - 1) + D^n^q^ + 



(150) 



By comparing the powers q^^q^ and q^ in this expression with the expan- 
sion of (149), we find, 

B|,| = -i (151) 

«m - ^^ (154) 

When a = 2, they reduce to the previous 13 = 1 result. These coefficients are 
written in terms of the inverse dimensional constants ^ ^ ,-,-, . The coefficients 

1 and a in (152) are the characters Xp(1) ^-s shown in the appendix A. 

34 



As we remarked, these coefficients -D[i,i], -D[a] and Dpi] are the coefficients 
of (146), and not of (147). For instance, the value of D[a] at /9 = 4 is different 
from the value of C[a] of (147) as 

The residual recursion equations, which gives the recursive relations between 
the coefficients, are derived for the coefficients of (147). However, we find 
that the following relations are also valid for the expressions of (151) - (154), 
which are the coefficients in (146), 

D[i] + (A; - l)D[i,i] + D[A] = (156) 

D[i] + (A; - 1)D[A] + 2D[n] = (157) 

The reason for the remarkable fact that the residual equations are satisfied, 
is that the residual equation is independent of the value of /9, and thus of the 
value of a. It means that the residual equation holds independently of the 
existence of the Vandermonde factor in (147). 

Therefore, for an arbitrary parameter a, we have found expressions for 
the coefficients which satisfy the residual equation. The coefficients C\gra-ph\ 
in (147) satisfies the same recursive equation. We have assumed that a = I. 
However, if we take this value of a as 



the coefficients Digraph] in (146) for the case of /5 become the coefficient C[graph] 
in (147) for the same /5. 

For instance, the previous D[a] becomes, after substitution of a = ^^, 

1 f-1 

which is indeed the value of C[a] in (147) as shown in Table B. Similarly, we 
find that -D[i,i] in (152) and Dpi] in (154) are the expressions given in Table 
B after the substitution of a = ^^, 

^'''^ i + a^k{k + a)^ k{k-iy k{{l-i)k-iy {l-i){k-iy 

(160) 
l + g (f-1) 

^["1 - 2kikTaj - A;((f-1)A:-1) ^^^^^ 
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The relation (158) is a duality relation, and by this relation, we obtain 
explicit expressions of the WKB expansion of the integral. This transfor- 
mation a = I — > 2Zfl becomes an identity for the case oi f3 = 1. For the 
value oi f3 = 2, this transformation becomes singular, and it needs special 
consideration (We will discuss this case in section nine). 

Ar order three, we have to consider the next order in (148). It becomes 

' (X[3,(l)^#«^ + X,,(l)^'^^'^'^^^^^'^'^ 



6(l + a)(l + 2a)^'"^^'^ ^ Z[3](I) '^'^'^' ' Z[2i](I) 

"l3](£)Z[l3](A) 



+X[i3](l)- 



Z,3](I) 



6k^^j2ik + a){l + 2af'^'^''^^^ 



"'^ -S3(£)S3(A) (162) 



3{k + a){k + 2a){k - l){k - 2) 



TT- .V. 1 f f~^ g(^ - g) f~, q\q-kf{2q-kf .^ ^ 

Usmg the values of S2[x) = ; , ssix) = — , it be- 

k k^ 



comes 

6k{k + a){k + 2a){k-l){k-2y ^ ^ ^ 

—Qakq^ + Qa^kq^ + SaPq'^ 

+12aq^ - 12a^q^ - 6akq^ 

+2g^ - 6ag^ + 2a^q'^ - ?,kq^ + 3aA;g^ + k'^q^) (163) 

When we put « = 2 in the above quantity, it becomes the same as (104). 
This quantity J*-^-* should be equal to the following sum of 8 terms, 

i^'^^ = -D[I,I,I] ■ [Ti2T34.T5q\ + D[Aj] ■ [ri2ri3r45] 

+ -D[Y] ■ [runsTu] + -D[II,I] • [T-i2r34] + -D[N] ■ [t"12T-13T34] 

+ ^[iii] ■ [ru] + ^[ Z ] ■ [^iVia] + ^[A] ■ [T-i2T-23n3] 

= -[D[i,i,i] ■ ^q'{q - lf{q - 2f + D[a,i] ■ q'{q - lf{q - 2) 

+ ^[Y] ■ lq\q - l)(g - 2) + D[n,i] ■ q\q - 1)' + Z^[n] ■ g'(g - 1)' 

+ D[in]-g' + /^[Z]-V(g-l) + AA]-0] (164) 

Comparing q^ term in (164) with (163), we have 

k^ + 3{a - l)k + 2{a^ - 3a + I) 
P'"''^ ~ k{k + a){k + 2a){k-l){k-2) ^ ^ 
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By the dimensional constants Zp{T) for general a, it is expressed as 

n, .- 1 r X[3](l) X[21](l) X[13](1) . . . 

where the characters X[p](1) ^-i^d the dimensional constant Zp{T) are given in 
Table D. (x[3](l) = l,X[2i](l) = -^— ,X[i«](l) = -^r- )• 

^ "i Cx ^ "i Cx 

By the dual transformation of (158), this expression becomes 

(^ _ 1)2 C! 

^PW] = -7770 77r\7770 7TI — ^[^ ~ 



A:((f-l)A:-l)((f-l)A;-2)^ {i - l){k - 1) 

which coincide, here also with the previous values listed in Table B. 

By the comparison of q^ terms in (164) and (163), we obtain the expres- 
sion of D[A I] as 

_ A: + 2a - 1 

^[^'^1 - ~A:(A; + a)(A; + 2a)(fc-l) ^^^^^ 

By the replacement of a = 2/(2 — P), it becomes 

f^-D^ 2 

^fAil = « ^ [1--5 1 (169) 

'^•^^ A:((fA:-l)((f-l)A;-2)^ (f-l)(A:-l)J ^ ^ 

which coincides with the result in Table B. 

By putting g = 1 in (163) and (164), we obtain the expression for -Dpii] 

as 

;i + a)(l + 2a) 

6A;(A; + a)(A; + 2a 
and by the dual transformation of (158), it becomes 

(^_1)(^_ 1) 
^™ = -*((f-l)*-l)((f -!)<:- 2) ^"'^ 

which agrees with the value in Table B. Note that for /9 = 4, this quantity is 
vanishing, and it means that for [3 = 4, there are no multiple line graphs. 

Extracting the contributions of _D[iii],D[a,i] and -D[iii] from (163) and 
(164), we have 

^[Y] ■ lq\q - l)(g - 2) + D[n,i] ■ q\q - if + /^[n] ■ q\q - if 
+D^ly2q\q-l) 



Dm = - X.^XJ::\ (170) 



q\q-l] 



■[3akq + llkq + 6q; g 



6A;(A; + a)(A; + 2a)(A;-l) 
+9ag-llg + 3aA;-7A;-6a2-15a + 7] (172) 
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Dividing both sides by g^(g — 1) factor, we obtain by putting q = 1, 

From the simple structure of [Y] (no multiple lines), we assume that 

1 ^ (f-1 ^^ 

k{k + a){k + 2a) ~ k{{^k - l)((f -l)k-2 



D,y. = i = ^-^-^ (174) 



which is also consistent with the known value for /3 = 4, /c = 4 (-D[y] = ^^ 
[4]). Then, we get 

1 + a 

^^^^^~2k{k + a){k + 2a) ^^^^' 

The sum of the remaining two terms is 

^W+^'-l = - M<. + a)(t!2a.)(^-l) l^'*-^"^ + °'+"P + °>'<™' 

Since -D[ii,i] is a coefficient of the term which has a double line in the graphic 
representation, it should be proportional to the factor (1 + a)/2, which cor- 
responds to the double line multiple factor. Therefore, the sum of (176) is 
divided into two parts, 

_ 1 + a , 2a , 

-^[11,1] - ~ nL.n. I ^\n. I n^\ l^ 



2k{k + a){k + 2a)' k-V 
1 .1 + a a(l + a). 



l + 2a^2Z[3](I) Z[2i](I) 



Dm = -777^— 4^^T7t[1+ " 



(177) 



k{k + a){k + 2a) k-1 

' ^' + "+.r^l (m 



(l + 2a)^Z[3](I) Z[2i]iiy 

These two expressions coincide with the values of Table B by the duality 
transformation (158). Since there is a cubic identity equation for three terms 
[ri^2'r34], [ti2T23'T34] and [ri2r23Ti3], the coefficients of these terms D[iij,D[n] 
and D[A] have ambiguities. We have to fix these ambiguities. Here we used 
a fixing assumption for the form of -D[n,i], namely that it has ^-^ as overall 
factor, and a pole at k=l. 

The last coefficient D[a] is not determined by the present choice of X and 
A. One method to determine D[a] is to use the residual equation, 

DA + ik- 2)D[N] + /^[A] + D[u,i] + /^[Z] = (179) 
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Using the previous expressions, we obtain 

1 a2 

^[^1 = ~k{k + a){k + 2a)^^ " r^^ ^^^°^ 

This reads to 

1 (1 + a)^ a^ 

This expression coincides with the value in Table B after the duality trans- 
formation (158). 

In Appendix D, we have derived the expressions of the coefficients D up 
to order fourth, from the extended zonal polynomial expansions. The results 
are shown in Table B. 



8 Series expansion for (3 = 4 



We have found in the previous section that the integral (1) may be expressed 
in two different dual ways, 

and 



h = H 



e-iE.<.-^. V - i Vy (i)^A^}M^] (183) 



I k 2 13 

where Xi = Xi — — ^a;^, and a = — , a' = — 1/(— — 1). The polynomials 
k -^^ p 2 

Zp{x) are the extended zonal polynomials with parameters a in (182) and a' 
in (183), and Xp(l) ^-^e the characters. 

In the previous section, the zonal polynomials were further expressed as 
polynomials in Tij as (146) and (147). 

The duality means that (182) and (183) give the same expression for the 
integral Ij3 under the relation of 
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or equivalently, 

/9 = - = 2-4. (185) 

a a' 

In the case /9 = 4, 6, ..., 2m for the even integers, a' becomes — 1, — |, ..., 
— ;^^. The zonal polynomial for a' = —1 was discussed before in a different 
context [16]. In such cases, the value of a' is negative, and the factor in the 

m— 1 

denominator ]^ (1 + qcn') in (183) vanishes, although the whole expression 

remains finite. Therefore, we need a special treatment for such (3 = 2m case, 
and we briefly discuss the case j3 = 4 here. The dimensional constants Zp{l) 
become degenerate in such cases. 

In the case (3 = 4, the parameter a' is —1. The dimensional constants 
show the degeneracy ; for instance, 2'[3](I) = k(k — l)(k — 2), and Z[i3](I) = 
k(k — l)(k — 2). Such degeneracies can be seen in the Table A. In addition 
to this degeneracy, the multiple line factors are proportional to (1 + a), and 
there are no multiple line graphs in /? = 4 in the representation of (183), and 
the expansion becomes rather simple for /5 = 4 as shown in [4, 5]. 

Here we study the zonal polynomial expansion in the form (183) for the 
(3 = 4 case, without writing it as a function of r. 

The expression (183) may be converted to more useful series in the case 
a = —1. In the previous section, we have written it in terms of the symmetric 
functions s„ for general a. Although there is a divergent factor jr^ in (183), 
this divergence is cancelled by the sum of the partitions of p. We investigate 
here the case a = — 1(/5 = 4) in higher orders, and obtain a useful expression 
for the integral /. 

We define (f){x, A) as 



perm. \.\.i<j[\-^i ^jjl^i ■^j)\ 



m=0 



m 



!n™V(i + 9«) p ' ^p(i) 



Note that we have f = e'' ^"^'^0(0;, A). Using the table of extended zonal 
polynomials and characters given in the appendix, we write the zonal poly- 
nomials in terms of the symmetric functions s„, and find the expression for 
0, by noting that ,51(5;) = 0. 



a _ ~ 
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+ 3(A; + a){k + 2a){k -l){k- 2) ^^W^s^A) 



+ 



a^ 



k{k - l){k - 2){k - 3){k + a - l){k + a){k + 2a){k + 3a) 

X [2ak'^{k'^ + ak — k + a)s4{x)s4^{y) 

+2ak{2k'^ + 3ak — 3k — 3a){s4{x)sl{y) + sl{x)si{y)) 
+ {k^ + 5ak^ - 5k^ + 6a^P - 18ak^ + 6k'^ - ISa'^k 
+18ak + 18a'^)sl{x)sl{y)] 

+ 0{x^) (188) 

There is no divergence in the hmit a ^ — 1 in this expression. Thus we 
obtain for /9 = 4, from the table in the appendix, up to order six , 

1 ~ k 

= 1 r, z7:S2(x)s2(X) H r. ZTrr, r:TS3(£)s3(A) 

^ 2{k-iy ^ ' ^ ' 3{k-lf{k-2Y -^^ > ^^ > 

[-2k'^{k'^ - 2k - l)si{x)si{\) 



8k{k-lf{k~2f{k-3f 

+2k{2k'^ -6k + 3){si{i)slC\) + sl{x)s4{\) 

+{k^ - lOP + 30k^ - 36k + 18)s^(£)s^(A)] 

+ 5(k^^Tnk^^Wik^^Wik^^ 

[S2{X)S3{X)S5{X) + S5(£)S2(A)S3(A)] 



(A;-l)2(A;-2)3(A;-3)2(A;-4)2' 

A;4 _ i4A;3 + 48A;2 _ 48A; + 24 ,., ,., ~, ~, 
+ 6{k - mk - 2Y{k - 3)2(A: - 4)2 ^2(a:)s3(x)s,(A).3(A) 

1 

^ (-5 + kf (-4 + kf (-3 + kf (-2 + A;)' (-1 + kf 

X {S6(£)S6(A)(--A;)[8 - 44 A; - 135 k'^ + 76 A;-'^ + 6 A;^ - 8 A;^ + A;^] 
6 

+ {s2{x)si{x)sQ{\) + S2(A)s4(A)s6(5))[-20 + 70 A; - 61 A;2 - 29 k^ 
+38k^ - II k^ + k"^] 

+ (S2(£)S4(£)S2(A)S4(A)) — [-2400 + 3600 k - 1740 A;^ - 240 k^ + 1099 A;^ 

8A; 

-708A;^ + 196A;6-24A;^ + A;^] 

+ (S3(5)^S6(A) + S3(A)^S6(£))77f80 - 200 A; + 251 A;^ - 284 A;^ 

6 

+154A;^-36A;5 + 3A;'^] 

+ (S3(5)'S2(A)S4(A) + S3(A)'s2(£)s4(£))(-777)[-400 + 200 A; + 525 ^ 

2k 

-690 k^ + 322 A;^ - 66 A;^ + 5 A;^] 
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+ (s3(£)^S3(A)2)^[-2400 - 1200 k + 7890 k'^ 
18k 

-8310 k^ + 4461 k^ - 1416 P + 264 k^ - 26 A;^ + k^] 
+ (S2(£)^S6(A) + S2(A)^S6(£))7^[240 - 769 A; + 882 /c=^ - 424 k^ + 90k^-7 k^] 

+ (S2(£)^S2(A)S4(A) + S2(£)^S2(A)s4(A))^[4800 - 11180 k + 11480 P 

-6775 k^ + 2384 A;^ - 481 A;^ + 50 A;^ - 2 A;^] 

+ (s2(£)3s3(A) + S2(A)'s3(£))t^[-1200 + 2270 k 

6k 

-1725 A;2 + 640 A;^ - 115 A;^ + 8 A;^] 

+ S2(5)^S2(A)^ ■ -i-[-25200 + 60960A; - 64030A;2 + 38192A;3 
48A; 

-13976A;^ + 3170A;^ - 431A;^ + 32A;^ - A;^]} 
+ 0{x'^) (189) 

We first check the simple A; = 3 case, given in the introduction. Muhi- 
plyiong by e~3(^i2+T23+ri3) ^-^ ^^ ^g recover the known result, 

^.l(^r,,+r,,+r,s)^ = 1 - ^ (^2 + r23 + ^3) 

1, , 

+ 7IT12T23 + T23T13 + T13T12J 

6 

- —n2r23Ti3 (190) 

where the last term is obtained from the identity 

-^(t-12 + T-23 + Tisf + ^(^2 + ^23 + ri3)s2(£)s2(A) 
1 ^ ~ 1 

+ -S3(a;)s3(A) = -— ri2r23ri3. (191) 

The fourth order term in (189) has a divergent coefficient —. However 

(A; — 3)^ 

the fourth order term is finite in the limit A; — > 3, and with the exponential 
term, it is vanishing at the end. In the fifth order, the same situation occurs. 
Therefore, we have recovered exactly (190) to all orders. 

The expression (189) is lengthy, however its large k limit becomes simple. 
In the large k limit, for each order x' , the coefficients of the products of the 

symmetric functions are of order -j. We take these leading terms, 

rv 

~ 1 - ^S2(£)S2(A) + — S3(5)S3(A) 
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1 ^ ~ I ^ ^ ~ ~ 

+ TJ^S5{X)S5{\) + —^S2{X)S3{X)S2{\)S3{\) 



48A;6 
+ 0(a;^) (192) 

There is a rule for the coefficients in (192). Only the combination of the 
same symmetric functions for x and A give the leading terms in the large 

1 
k limit. For instance, S2{x)'^S2{X) and S2{x) S4{x) S2{X) S4{X) are of order — . 

The coefficient of S2(x)'^S3(A)^, which has different symmetric functions, is of 

order — -. 

The coefficients of (317) is obtained as follows. For the term of s„(a;)^Sm(a;)*, 
{n 7^ m), the coefficient becomes 

For the general case, s^\s'^l ■ ■ ■ Snj, the coefficient is proportional to 

. The order of x is given by n^^ + ■ ■ ■ + n^^ , and if 



Pi\p2^- ■ ■ ■pj\n\^n!2 ' ' '^^^ 

this order is even, the minus sign has to be included. Thus, we have the large 

k expression for 0. 

In the large k limit for fixed Xj, A^, / is given by e^ fc ^ ^'^ , and 0(x, A) = 1, 
for all values of a, as shown in Appendix D. However, in some problems, in 
which Sn{x) is not order of one, the above large k formula might be important. 

9 The HIZ-integral for (3 = 2 and the char- 
acter expansion 

We use here the same zonal polynomial method as for (3 = 1. We define 
^a = Xa — \ Z]b=i ^b- The integral becomes (a = | = 1) from the expression 
(146), 
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This series expansion has to reduce to the simple closed form 

b^2= E -TT—^ (195) 

perm.ofXi '■'■i<j *J 

in order to reproduce the original HIZ formula [11, 12].. 
In other words one must prove the identity, 

^ eJE^. ^ ^ ^ 1 ^1 Z,{x)Z,i\) ^^gg) 

perm. H^ij m p "^' 11^0 (^ + ?) Zp{l) 

The proof of this identity is easily done by writing the zonal polynomials, 
which are Shur functions, as the ratio of the determinants, 

Zp{x) = ^^ (197) 

det [xj ] 

The product of the determinants is also a determinant, it is simply det[e^''^j], 
from the Binet-Cauchy theorem [18]. 

This identity leads to interesting equations. The right hand side of (196) 
is written by the terms of r, as shown in the case oi f3 = 1. 

For the case k=2, the proposition (196) is easily proved by expanding the 
exponent. For k=3, we have the following identity : for instance, 

E ^"'' ^ "'' ^ "''^' = (198) 

perm. T12T23T13 

for p=l,2 and 4. Let us remind the reader here, that the permutations in 
this sum are interchanges of the Aj's for fixed Xj. For p=3, we have 

E ^^'' ^ ""'' ^ ^''^' = 81 (199) 

perm. T12T23T13 

For p=5, it becomes 

perm. Tl2T23^13 4 

= 405 [(rf2 + r^3 + r^s) - {ri2r23 + r^ns + T23T13)] 

(200) 

where the second equality comes from the expression of S2{x) of (95). For 
p=6, we have 

E ^"" ^ "'' ^ "''^' = iSl)%{x)s3{X) (201) 

perm. T-i2r23ri3 
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For p=7, 

E ^"" "" "'' '' "^^^' = 1701(^)^h(£)s,(A)]^ (202) 

perm. T12T23T13 ^ 

where for k=3, we have 54(0;) = |[s2(5)]^. For p=8, 

E ^^'" ^ ^'' ^ '"''^' = 648 X &'[s,{x)ss{x)s,C\)ss{\)r- (203) 

perm. ^12T"23l"l3 ^ 

From these resuhs, one checks the identity (196). 

To understand these identities, we divide the left hand side of (200) into 
five different types of terms. They are all represented by S2(5)s2(^)- 



perm. 



Tl2T23n3 4 



^ rf^ins + r23) + ■ ■ ■ 225 ~ 

perm. T-12T"23T'l3 4 

E (^12 + 4 + 4) = 18S2(£)S2(A) 
perm. 

E(''"l2'''23 + ■ ■ ■) 9 /~\ /Tx 
— — = -;S2{X)S2{X) 



perm. 



T'12T"23T"13 4 



E (^12^23 + ■■ ■) = IS2{X)S2{\) (204) 



perm 



After summing over permutations, the r-series are expressible by symmetric 
functions of x and A, and are given by homogeneous polynomials in r. 

We thus have shown that the zonal (character) polynomial series can give 
the expression for /, although / is just one. 

The dual representation of (8) is singular for f3 = 2. The transformation 
a = 2ZI3 becomes divergent at /3 = 2. However, if we write / as a r expansion 
for fixed k, and let a — * 00, we have an infinite series, as can be seen in Table 
B. There are non-vanishing terms in the limit a = 00. For instance, at first 
order, Ci = —^. This looks quite strange, since we know that / is one. In 
the following, we consider this puzzling problem. 

Collecting the no n- vanishing terms in the limit a ^ 00 from Table B, we 
find the following expression, 

perm. of Xi -l -1 *J V / 

+ j^{r'r2 + ---) + 0{x% (205) 
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According to this expression, /^=2 is not one. This discrepancy may be 
understood by looking at the simple example k=2. In this k=2 case, we have 

h=2 = E (l-^n2 + Tri'2 + ---) 

perm.ojAi ^^ 

^a::iAi+3;2A2 ]^ J^ ga::iA2+a;2Ai j^ j^ 



(1 - ^5^2 + Tri^2 + ■ ■ ■) (1 + 15^12 + -^rt, + 

Ti2 2 4 ri2 2 4 



(206) 



This expansion is divided into two parts; 





ga;iAi+a;2A2 _ ga;iA2+a;2Ai 


/:l=2 


ri2 




- ^[e^^'^+^^'^(l-^ri2H 



■) + e"^^^+"^^Hl + iri2 + ■ ■ ■)](207) 

Remarkably, if one expands this expression in powers of the x's and the A's, 
one finds cancellations between the first and second term between brackets, 
and agreement with (195). 

Thus we find that the limit f3 = 2 is somehow singular, and the r- 
expansion, with the coefficients of Table B for f3 = 2, gives one for the 
value of the function /^=2 after summing over the permutations. 

We have thus recovered the ordinary HIZ formula of (195). The r- 
expansion for (3 = 2 oi (205) may be useful for a check of the coefficients 
for general (3, since it should reduce to one for fp=2- 

10 Expansion of fp for large (3, and for large 
k and fixed (3 

We note that the final result for ff^ in the large f3 limit is simple. The 
coefficients C involve the products of the inverse of the multiplicity /!. For 
instance, the term of T'^2'^24'^23'ri3T34 has a coefficient as 

^1^2^23^-13^-34 (208) 



2!3!A;8 
Namely fp is 

/. = E|;^^(-i)"n-.,, (209) 

where n is the total number of the bonds Tij, and Ij is the multiplicity. The 
expression (209) is simply equl to 

f^ = e-^^^<j^'^ (210) 

46 



This formula is valid for large f3 and for arbitrary k. It may be interesting 
to refine above expression. For fixed f3 and in the large k limit, we have an 
expansion in powers of r. We may thus expand (8) for large k. / is given by 

and as shown in Appendix D, or shown in the Table B, it is transformed into 
a r-series. We find that jp is given by 

Of course it reduces again to (210) in the large /9 limit. The above expression 
(212) automatically satisfies the condition that the series of fp stops at the 
order | — 1 in r when /5 is an even integer. Thus we have obtained the 
improved expression of (212) for /^, which is valid in the large /5 limit or in 
the large k limit. This asymptotic form (212) may be useful for finite fixed 
/? (for instance, /9 = 1) and large k. 

In the Appendix D, the large k limit of the coefficients C at /-th order is 
investigated. They are given by 



^' I (l+0(-^ 



C = (-i)'^Ti^rT^— t(i + OO) (213) 



where 5" is a degeneracy factor for the multiple lines. It is remarkable that 
the only dimensional factor which appears in the large k limit, is the first 
row of the Young tableau, which has the form of {k + ma). For instance 
(/c — 1),(A; + q; — 1),.., which are the second row factors, do not appear in the 
large k limit. 

11 Summary and discussions 

In this article, we have given the expression of the integral (1) as series in 
the variables Tij = {xi — Xj){\i — \j) for the function / defined in (2). As 
discussed in section seven, from the extended zonal polynomial expansion 
of the integral J, and the use of the dual representation of (158), we have 
obtained expressions for this same function /. The coefficients Digraph] of 
this r expansion are expressed through the dimensional constant Z[p] (I) of the 
extended zonal polynomials (Jack polynomials). The results for these C[graph] 
coincide with the direct perturbational calculations developed in the section 
five, where the residual equations among various coefficients C[graph] are used. 
It is remarkable that these recursive residual equations are indpendent of the 
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parameter f3. We have found the exphcit expression for the WKB expansion 
by this duahty. 

In this paper, we have proved that the extended Harish-Chandra-Itzykson- 
Zuber integral for the general (3 case is expressed by the variables of tij. The 
proof is given in two stages; the first is the expression of the Harish-Chandra- 
Itzykson-Zuber integral by the zonal polynomial expansion with the param- 
eter a, which is i^^ by the duality. The second is the transformations of the 
products of the symmetric functions s„(x) in the zonal polynomial expansion 
into the r variables. This transformation is discussed in Appendix B, and in 
Appendix D in detail. We found that there are identities among the r terms. 
For instance, we have found explicitly these cubic and quartic identities in 
the appendix C. These identities give a sort of gauge freedom to choose the 
values of the coefficients of the r expansion. We have considered a fixing 
of these ambiguities from the large k behavior by imposing definite asymp- 
totic forms. For the /? = 4 case, only single line graphs appear, and for this 
reason, there are no ambiguities for the coefficients C, which are uniquely 
determined. 

The integral (1) is important for the investigation of the random matrix 
theory, specially in the presence of an external matrix source, as shown for 
the 13 = 2 m. [3, 7, 8, 9]. We will discuss in a separate paper the applications 
of the present results [10]. 
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Appendix A: Extended zonal polynomials (Jack polynomials) 

We have uses the expansions of the extended zonal polynomials ( named 
Jack polynomials) in section seven, and we have derived the expression for 
/ from the extended zonal polynomial expansion of the HCIZ-integral Ip by 
duality. Therefore, in this appendix, we give the needed important quantities, 
the characters Xp(1)) ^"^^ the dimensional constants Zp{\) [13, 14, 16]. 

The lower Jack symmetric polynomials Z\p^{X) and their dimensions 



Z[p 


l(I) 


are 






















%]W = 


-- Si, Z[i 


|(I)=k, X[i](l) = 


1 














Z[2]{X) = 


-- si + as2, 


^[2] (I) 


= k(k + 


«), 


X[2] 


(1) = 


= 1 








ZiHX)- 


= S\- S2, 


%^](I) 


= k(k- 


1), 


X[l2 


](!) = 


= X[i2 


](1) 






%]W = 


-- s\ + 3asiS2 


+ 2a2s3, 















a 



48 



Z[3](I)=k(k + a)(k + 2a), X[3](l) 



Z[2i]{X) = .sl + {a- l)siS2 - as3, 



Z[2i](I)=k(k + a)(k-l), x[2i](l) 



Qa(l + a) 



2 + a 



Zm3i(X) = sl- 3siS2 + 2s3, ^[13] (I) = k(k - l)(k - 2) 



X[i3](l) 



a2(l + 2a) 
2 + a 



(214) 



where a = [3/2. The classical symmetric functions are denoted by s„, s„ 
Then one has the relation, 



(trX)^ 



nm=\(l + ma) 



[EXp(l)Zp(x)] 



(215) 



where p is a partition of the integer g, when a = 2/(3 = 2. 

The dimensional constants Zp{T) are obtained by putting X = / in the 
zonal polynomials Zp{X). They are factorized as a polynomial in k. 

From the constants given in (214), and the sum rule (215), we find the 
HIZ-integral for general j3 : 



T - ^ 1 1 V,. .,. Zp{X)Zp{A) 

' i'o-!n-V(i + g«)t^^^'^ z,(i) 



(216) 



If A=I, it becomes 

Ip = e*""^ (217) 

which is the correct expression by definition. The values in the following 
tables of coefficients of the zonal polynomials come from [14] (with a minor 
correction). The characters Xp{^)j ^-^e then evaluated on the basis of these 
values. It agrees with [15] when a = 2 up to sixth order. 



The extended zonal polynomial (Jack polynomial) Zp(a;) with a 
parameter a and its coefficient of the symmetric functions 

Xp(l) is a character and Zp{l) is a dimensional constant. 

/ = 1 





Si 


Xp(l) 


m) 


Z[i] 


1 


1 


k 



49 



1 = 2 












si 


S2 


Xp(l) 


z,{i) 


Z[2] 


1 
1 


a 
-1 


1 
a 


k{k + a) 
k\k - 1) 





si 


S1S2 


S3 


Xp(l) 


Zp(S) 


Z[z\ 

Z[2,l\ 
Z[13] 


1 
1 
1 


3a 

a-l 

-3 


2a^ 
—a 

2 


1 

6a(l+a) 

2+a 
a^{l+2a) 

2+a 


k{k + a){k + 2a) 
k{k + a){k-l) 
k{k-l){k-2) 



1 = 4 
















sf 


S?S2 


si 


S1S3 


S4 


Xp(l) 


Z[4] 




6 a 


3a' 


8a2 


Ga^' 


1 


Z[3,l] 




3a -1 


—a 


2a2 - 2a 


-2a2 


6a(l+2a) 
1+a 


^[22] 




2a -2 


a^ + a + 1 


-4a 


a — a^ 


6a^(l+3a) 
{l+a){2+a) 


Z[21^] 




a-3 


—a 


2 -2a 


2a 


6a'^(l+2a){l+3a) 
{l+a){3+a) 


y 




-6 


3 


8 


-6 


a-^{l+2a)(l+3a) 


Z/[l4] 




(2+a){3+a) 





z,{\) 


Z[A] 


k{k + a){k + 2a){k + 3a) 


Z[3,l] 


k{k + a){k + 2a){k-l) 


^[22] 


k{k + a)lk + a - l){k - 1) 


Z[21-2] 


A;(A; + a)(A;-l)(A;-2) 


Z[14] 


k{k - l){k - 2){k - 3) 



1 = 5 





s\ 


S?S2 


sisl 


sfsa 


S2S3 


S1S4 


S5 


Z[5] 




10 a 


15a2 


20a2 


20a3 


30 a^ 


24 a^ 




Z[4,l] 




6a -1 


3a(a-l) 


a(8a - 3) 


-W 


6a2(a-l) 


-6a3 




Z[32] 




2(2a - 1) 


3a2 - a + 1 


2a(a - 3) 


2a{a' + 1) 


-a(7a-l) 


-2a2(a- 


1) 


Z[31'2] 




3(a-l) 


-5a 


2{a-lf 


-2a(a - 1) 


-4a(a - 1) 


4a2 




Z[22,l] 




2(a-2) 


a^ - a + 3 


-2(3a-l) 


-2{a^ + 1) 


-a{a - 7) 


2a(a-l) 




Z[2l3] 




a-6 


-3(a-l) 


-(3a -8) 


5a 


6(a-l) 


-6a 




Z\i^] 




-10 


15 


20 


-20 


-30 


24 
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Xp(l) 


^p(I) 


^[32] 
^[312] 
^[22,1] 
^[213] 


1 

20a(l+3a) 

2+3a 
30a2(l+4a) 

(l+a){2+Q) 
30a2(l+2a){l+3a){l+4a) 

(l+a)(3+2Q)(2+3a) 
30a2(l+3Q,)(l+3Q,)(l+4a) 

(l+a)(2+Q)(3+a) 
20a3(l+2a){l+3a){l+4a) 

(2+o)(4+a){3+2a) 
a'^(l+2a){l+3a){l+4a) 

(2+a)(3+a)(4+a) 


A;(A; + a)(A; + 2a)(A; + 3a)(A; + 4a) 
k{k + a){k + 2a){k + 3a){k - 1) 
k{k + a){k + 2a)(A; - l){k + a - 1) 
A;(A; + a)(A; + 2a)(A;-l)(A;-2) 
k{k + a){k - l){k + a - l){k - 2) 
A;(A; + a)(A;-l)(A;-2)(A;-3) 
A;(A;-l)(A;-2)(A;-3)(A;-4) 



1=6 





4 


S1S2 


„2 2 
■^1-^2 


4 


sfss 


Z[6] 




15 a 


45a2 


I5a^ 


40a' 


Z[5,l] 




10a -1 


3a(5a - 2) 


-3a2 


4a(5a - 1) 


Z[42] 




7a -2 


9a^ - 5a + 1 


a(3a^ + a+ I) 


8a(a- 1) 


Z[32] 




3(2a-l) 


3(3a2-a+l) 


-(5a2 + 3a + l) 


4a(a - 3) 


Z[4l2] 




6a -3 


3a(a-4) 


-3a' 


2(4a2 - 3a + 1) 


^[321] 




4(a-l) 


3(a-l)2 


—a{a — 1) 


2a' - 9a + 2 


2'[313] 




3a -6 


-3(4a - 1) 


3a 


2(a2-3a + 4) 


^[23] 




3(a-2) 


3(a2-a + 3) 


a(a2 + 3a + 5) 


-4(3a - 1) 


Z[2212] 




2a -7 


a^ - 5a + 9 


-(a2 + a + 3) 


-8(a - 1) 


-Z'[214] 




a -10 


-3(2a - 5) 


3a 


-4(a - 5) 


^fl«l 




-15 


45 


-15 


40 



(continued) 





S1S2S3 


4 


S1S4 


S2S4 


Z[6] 


I20a^ 


40 a^ 


90a3 


90a^ 


Z[5,l] 


20a'{a-l) 


-8a3 


6a2(5a - 2) 


-18a3 


Z[42] 


4a(2a- l)(a-l) 


-2a2(a - 1) 


a(6a2- 17a + 1) 


a'{6a'-a + 5) 


^[32] 


12a(a2 + l) 


2a'{2a' + 3a + 3) 


-3a(7a - 1) 


-3a(4a2 + a + 1) 


Z[Al2] 


-6a(3a - 1) 


4a' 


6a(a-l)2 


-Qa'{a - 1) 


Z[321] 


(a-l)(a-2)(2a-l) 


-a{2a' + a + 2) 


-9a(a-l) 


-2a(a-l)2 


Z[31-i] 


-6a(a - 3) 


4a' 


-6(a-l)2 


6a(a-l) 


^[23] 


-12(a2 + l) 


2{3a' + 3a + 2) 


-3a(a - 7) 


-3a(a2 + a + 4) 


Z[22l2] 


-4(a-l)(a-2) 


2a(a-l) 


-(a2-17a + 6) 


5a2 - a + 6 


^[214] 


20(a - 1) 


-8a 


6(2a - 5) 


-18a 


Z[16] 


-120 


40 


-90 


90 



(continued) 
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S1S5 


■% 


Xp(l) 


Z[&] 


144a4 


120a^ 


1 


%,1] 


24a3(a_l) 


-24a^ 


15a(l+4a) 
l+2a 


^[42] 


-4a2(5a-l) 


-Qa^{a - 1) 


90a2(l+2a){l+5a) 
(l+a)2{2+3a) 


^[32] 


-12a2(a - 1) 


-2a'^{a-l){2a- 1) 


30a3(l+4Q)(l+5a) 
(l+a)2{2+Q)(l+2Q) 


Z[4l2] 


-12a2(a _ 1) 


12^3 


10a2(l+3Q,)(l+4a)(l+5a) 
(l+Q)2(l+2a) 


^[321] 


-a(2«2_i3^^2) 


4a2(a _ 1) 


720a-^(l+a)(l+3a)(l+4a)(l+5a) 
(2+a)2 {l+2a){3+2a){2+3a) 


%13] 


12a(a-l) 


-12a2 


10a^{l+2a){l+3a)(l+4a){l+5a) 
(l+a)2(2+Q)2 


^[23] 


12a(a-l) 


2a(a-2)(a-l) 


30a4(l+3Q,)(l+4a)(l+5a) 
(l+a)2{2+Q)2(3+a) 


Z[22l2] 


4a(a - 5) 


-6a(a-l) 


90a4(l+2a)(l+3a){l+4a){l+5a) 

(l+Q)2(3+a)(4+a)(3+2Q) 

15a4{l+2a){l+3a){l+4a){l+5a) 

(2+a)2(3+Q)(5+a) 


^[214] 


-24(a - 1) 


24a 


y 


144 


120 


a5(l+2a)(l+3Q)(l+4a)(l+5a) 


Z/[l6] 






(2+Q)(3+a)(4+a)(5+a) 












^p(I) 




^[6] 


A;(A; + a){k + 2a){k + 3a)(A; + Aa){k + 5a) 






%,1] 


fc(A; + a)(A; + 2a){k + 3a)(A; + 4a)(A; - 1) 






^[42] 


k{k + a)(A; + 2a){k + 3a)(A; - l){k + a - I) 






^[32] 


k{k + a){k + 2a){k - l){k + a- l){k + 2a 


-1) 




Z[412] 


k{k + a){k + 2a){k + 3a){k - l){k - 2) 






-Z'[321] 


klk + a){k + 2a){k - l){k + a- l){k - 2) 






-Z'plS] 


k{k + a){k + 2a){k - l){k - 2){k - 3) 






^[23] 


k{k + a)(A; - l){k + a- l){k - 2){k + a - 2) 




Z[22l2] 


k{k + a)(A; - l){k + a- l){k - 2)(A; - 3) 






^[214] 


k{k + a)(A; - 1)(A; - 2){k - 3)(A; - 4) 






^fl«l 


k{k - l){k - 2){k - 3)(A; - A){k - 5) 







Appendix B: The transformation of the paired products of sym- 
metric functions Sn{x) (power sum) to r-polynomials 

We consider the transformation of the paired products of the classical 
symmetric function s„(x), and s„(A) to the r-polynomials, where Xj = Xj — 



TL^Xj, and Tij 



[Xi 



Xj){\i 



\j), and Snix) 



E^r- 



At oprder two, S2{x)s2{X) is expressed in terms of t^j by 



S2{X)S2{\) 



^11 ■ F, 



tjl 



O-A ■ [runsl + ^1,1 ■ [^-12^34 



(218) 



where the coefficients a are functions of k. We apply the differential operators 



k,l dxidxjdXkdXi 



on both sides of the above equation. We obtain 

^2,2 [^12] = 4, dIiI[ti2Ti3] = £'2,2 [^12^^34] = 



(219) 
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We use the notation [rfg] for J2i<jTij, [T12T13] for J2i<j<k''~i,j''~i,k,Gtc. For the 
symmetric function S2{x)s2{X), we have 

Dliis2ix)s2i\)) = ^^^^ (220) 

Therefore, we obtain the coefficient of [1^2] as an = ~^2 ■ Similarly, we 
apply 1^2 '3 s-'^d -024) S'nd then we find the coefficients of [ri2ri3] and [ri2r34], 

DlilrU = 0, Dl:l[n2m] = 2, Dli[r',,] = 0, Dliis2ix)s2i\)) = "^^^^ 

(221) 
These equations give the coefficient of [ri2ri3] as a a = — ^5"^- 

^2,4 [^^2] = 0> ^2,4[n27"l3] = 0, /^2,4[n27-3,4] = 2, Dl'^^{s2ix)s2iX)) = J^ 

(222) 
These equations give the coefficient of [ri2r34] as o"ii = -p-. These results give 
(95). 

We have examined all possible differential operators of order two for (218) 
; they are D^, DH, DI2, Dl\, D^, DH, DH. These operators con&m the 
equation (218) with the coefficients determined hereabove. ' We have thus 
estabhshed the identity (218). 

At order three, we transform S3 (5) S3 (A) to r-polynomials. 

S3 (X) S3 (A) = CTni ■ [t^2] + ^L ■ [^iVls] + O-A ■ [n2T-13T"23] 

+ O^Y ■ \T\2TViT\A\ + f^N " V\2TViT'ii\ + 0-n,I " [t-iV34] 

+ O-A,! ■ [T-12n3^45] + 0-1,1,1 ■ [n2l"34^56] (223) 

By applying -02,2,2 = dSdx^ '-"^ both sides of (223), we find two non- vanishing 
contributions : 



k^ 

From this result, we obtain am = _i±:il_|±:^_ From -D2,'2,'3 5 W6 obtain 



Dlii[r'^2] = -36, /^S(S3(£)S3(A)) = 36^ ^^ '- (224) 

(fc-i)^(fc-2)2 p_^ ni.1'1 



^2;2;3KVi3] = -12, Dl:liis,ix)ssix)) = sq^-L^^^^ (225) 

which gives a^ = li" ■ By the differentiation -02,3,4 5 we obtain 

Dl'kllrurisru] = -6, DI:1:1{ss{x)ss{\)) = 72 ^^~^^f ~^^ (226) 
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which reads to cry = —^{k — l){k — 2). For the differentiation -D2,'3,'5! noting 
that [Ti2ri3r45] includes a sum of the relevant terms T12T1ZT45 + t^ti^tu + 
T13T15T24, we obtain 

DliilruTisU,] = -6, Dl:l:t{s,{x)s,{X)) = -72^ (227) 

which reads (Ta,i = 12^p^. By the differentiation D2\'q, we obtain 

Dl'Simrun^^] = -6, Dl'',ll{ss{x)ss{X)) = 144^ (228) 

which reads aijj = "p- Thus we determine 5 coefficients of a as the func- 
tion of k, uniquely. These coefficients are represented in (101). For other 
3 coefficients, o"a, o"n, o"ii,i, we need other differential operators, which give 
coupled equations. 

By the differentiation D2 2\ of (223), two terms of r are non- vanishing. 



D2\2,aW^ • [^"12^-13^24] + 0-11,1 ■ [^"12^34]} = -4o-N - 4aii,i 

r 



Dlii{ss{x)s,{X)) = 36^^-^ (229) 



which reads 

(k — 2)'^ 
ctn + ^ii,i = -9 ^ ^4 ' ■ (230) 

By the differentiation -03,3,2, we obtain from (223) 

(k — 2)'^ 
8a^ - 4aA + 4(A: - 3)an,i = 36 ^ ^/ (231) 

Using the value of a^ , we have from above equation, 

(15-6A:)(A;-2)^ 
~aA + [k- 3)crn,i = 75 • (232) 

By the differentiation -D2'2'3; ^^ obtain 

(T^ + a^ + {k- 3)aN = 9 ^^~, (233) 

Using the obtained value of a^ , we have from above equation, 

^(1. ',^ 3(A;- 2)^(4 -fe) 
aA + [k- 3)o-N = 7i • (234) 
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The three relations of (230), (232) and (234) are not hnearly independent. 
Indeed if we shift anj -^ anj + a, o"n — ^ (Tn — « and a a — > (Ta + ^(^ — 3), 
these relations are unchanged. It means that there is a cubic identity, 

/3 = [II,I]-[N] + (k-3)[A]=0. (235) 

The proof of this identity is discussed in the appendix C Therefore, the 
coefficients (Tn,i, ctn, o"a have ambiguities up to a parameter a, which is an 
arbitrary constant. One can add a term al^, i.e. zero, to the expression of 

S3(5)S3(A). 

This may be used to write simple expressions for these coefficients, for 
instance, 

(A; -2)2 (A; -2)2 (A; -2)2 ^ ^ 

aii,i = -6^^, ^N = -3^^, ^A = 3^r^ (236) 

We have examined the transformation of (223) by all possible differential 

operators Dlt'n on (223). These operators are, DlH, DlH DlH DlH DlH DlH 
nii3 ni23 nii2 nii2 n^ii n^u n^u ni23 nii2 ni23 ni23 ni23 nii2 rpi • 

-^223 ! -'^123 5 -'-^123 ) -^^223 ? -^234 5 -^334 5 -^234 ) -'^124 ! -'^134 5 -^145 ? -^456 5 -'-^124 ) -^^345 • -■- ^io 

means that we have proved (223) with explicit coefficients. 

In fourth order, the products of symmetric functions {s2{x)s2{\))^, 
S2(x)2s4(A) + 54(5)52(^)2, and S4(£)s4(A) appear in the zonal polynomial 
expansions. Since S2(5)s2(A) is transformed as (218), (s2(5)s2(A))2 is trans- 
formed as the square of (218). 

We first consider the transformation of S4(a;)s4(A). It is written as a sum 
of 23 possible terms, 

S4(£)s4(A) = 0-[+][ri2ri3ri4ri5] +(T[^^][ri2ri3r34r35] 

+ 0-[>] [t'12T-13T"14T"23] + (^[n] [t-12T-13^24T"34] 
+ Cr[M] [ri2T-i3r34r45] + 0-[n,I] [ri2ri3r34r56] 

+ 0-[Y,I] [n2n3n47"56] +0-[A,I][Tl2n3T"23T45] 
+ 0-[A,A] [t"12T-13^45''46] + (^[A,!,!] [^"12^13^45 ^6?] 
+ 0"[I,I,I,I] [t"12T34T56T"78] + ^^I U ] [t-iVi3T"24] 

+ o-Q] [r22ri3r34] + a^ a ] [t"iVi3T"23] 

+ o-[ Z,i] [^iVi3T45] + o-[A,ii] [rurnrl^] 
+ cr[n,i,i] [TuTune] + o-[ii,ii] [t'iV34] 

+ qiiii][4] (237) 

A systematic way to determine these coefficients a consists in classifying 
the terms by the number / of the points in the graph. 
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For / = 8, we make use of the differential operator D^jg, and the graph 
is [I, I, I, I]. Then, we find uniquely that 

216 , , 

^M,M = ^ (238) 

For 1=7, we have two differential operators, DlHj and -Disgy, which act on 
the possible two graphs [I, I, I, I] and [A, I, I] (other graphs have less than 7 
points and do not contribute). From -DI5I7, we obtain 

(TA,i,i = -72^^^ (239) 

For 1=6, in addition to the above two graphs, we have 4 graphs, [N, I],[Y, I], 
[A, A], and [II, I, I]. We have 6 different differential operators, Dllll,DllJl, 
DllllDllll, Dllll^nd Dim From DHH we find uniquely, 

^[Y,i] = 72 (240) 



From -D34II, we find 



From -DI4II, we have 



,(2^3) 
A;6 



^[A,A1 = -72^^\r^. (241) 



-24aii,i,i + 12(Ty,i - 48(Tn,i - 12(A; - Q)aA,i,i = -24 ■ 24 ■ S^-r^ (242) 

which reads to 

(2A;2 - lOA; + 15) , , 

^ii,i,i + 2aN,i = 36^ ^ (243) 

From -D2256! ^^ obtain the same relation. 
From -DI256' ^^ obtain 

8an,i,i+14(A;-6)aA,i,i+2(A;-6)(A;-7)ai,i,i,i+4aA,A+16aN,i = 24-24-9— (244) 

which gives again, 

(2A;2 - lOA; + 15) , , 

^11,1,1 + 2aN,i = 36^ '- (245) 

Thus we are unable to find definite valuesfor o"n,i,i and o"n,i. The reason of 
this indefiniteness is the existence of an identity, which is an extension of 
the cubic identity. In the appendix C, we have derived the following identity 

(282) : 

[II,I,I]-[N,I] + (k-5)[A,I] = 0. (246) 
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This identity gives a freedom in the choice of their respective values. 

At next order / = 5 (five points), in addition to these 6 graphs, we have 
to consider six new graphs, [X],[ZZ],[M],[A, I],[Z, I], [A, II]. 

For the differential operators D"-J}f^, we have the following different kinds 
for / = 5, 



7-)llll nlll2 7-)1122 7-)1123 nlll^ T-)1122 

-^2345 5 -^2345' -^2345' -^2345' -^3345' -^3345 

r-)1234 r-)1123 f)U23 nll24 r)^l34 nH^^ 

-^1235' -^1145' -^1245' -^^2235' -^2235' -^^1345 



(247) 



where we note that DHH = DHH = DHHMEI = D\fi = D\fi. (the 
equal sign means the equivalence for the operator in this problem). 
We have uniquely from -D2345) 

(A;-l)(A;2-3A; + 3) , , 

^x = -72^ ^^-p '- (248) 

We next obtain the coupled equations. From -D33II, we have 

(A;-3)(2A;-3) , , 

(TM + (TA,u = 36^ ^^ '- (249) 

From -D2345! ^^ obtain 

ax+{k- 5)(TY,i = -288 ^^'"^^^ + ^^ (250) 

rv 

From -D334I, we obtain 

(A;-3)(A;2-3A; + 3) , , 

^11 + 2a^, = -48^ ^^^ ^ (251) 

From Dim, ^^ have 

ax - Q(Tii - 3(A; - 5)aY,i - 12a^ j = 288^ ^ (252) 



which reduces to (251). 
From -D2345! ^^ have 



2(JA,n + 2cTM + (A; - 5)aA,A = 144^^^—^ (253) 



This reduces to (249). 
From Dllll, we get 



A 



^11 - 60-M - 2(A; - 5)crN,i - 2o-A,i - {k - 5)aA, 

-2a^ J - 2aA,ii = 144 ^^"/^ (254) 
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From -D22355 ^^ have 

(k — 3)^ 
-au + 2o-^^ + {k- 5)o-N,i - cta,: + 2(T/ j +{k- 5)o-ii,i,i + aA,ii = -144 — — — 

(255) 
From -Djilf, we have 

ax - 8aii + 80-M + 8(A; - 5)o-n,i - 4o-y,i + 4o-a,i 

+ (A; - 5)(A; - 6)aA,i,i - Sa^j + 2(A; - 5)an,i,i = 144 ^^'/^ (256) 

These three equations coincide when we put the known values of the 
coefficients. 

The last differentiation in (247), -D1235' gives 

ax - ^'^crii - l2aM - 2A{k - 5)o-n,i - 4(A; - 5)o-y,i 

-9{k - 5)cta,a - 12(A; - 5) (A; - 6)aA,i,i - (A; - 5) (A; - 6) (A; - 7)cti,i,i,i 

-12aA,ii - 12(A; - 5)aii,i,i - 24ct^ j = 5184-^ (257) 

However, this leads to known identities rather than to a new relation. 

Several coefficients of the terms, which have less than 4 points, are deter- 
mined uniquely by the differential operators. We list them here, 

qmi] = l-,{k-lf{e-3k + 3f, {dI'HI) 

^k-l)ik-3){e-3k + 3) 1,1,1,1 

^[H = ^2 — , (i^2,2,3,4J 

(fc-i)(2A:-3)(A:^-3A: + 3) 1,1,1,1 

^[«] = ^ p ' 1^2,2,3,3) 



(A;-l)(A;2-3A; + 3) 



2 



-[Zl = -4^ ^li -^ (D 



A;6 



'^2,2,2,3J 



(2 A; - 3) (A;^- 3 A; + 3) 1113 

(Jd = -12 — , (1^2244)- (258) 

We have coupled equations by other differential operators. We list them 
with the differential operator Df, which are used. 

apn,!] + a[ n ] = 16 ^^''^^^ + ^^' , (D^,4'ii) (259) 

-2(7 [^] - (A; - 3)a[iii,i] + a[A] = 16 -^ (^2,2,2,3) • (260) 

(/^ _ ^)(2k — 3) 
3a^ + 2a[u,u] + ap] + {k - A)aM + (k - 4)(7a,ii = -72^ 



A;6 

2,2,3,4J 



[DliiD- (261) 
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,(2A;-3) 



2 



2a[n,n] + ap, = 36^—-^, {D'^i%). (262) 



-2(A; - 3)a> + 6a« - Qa^ + {k - 3){k - 4)aA,ii + 2(A; - 3)aii,ii 



(2A;-3) 



2 



+4(A; - 3)(73 = 72^-^, (Z^^^a.ts)- (263) 

ik-A)a,, + a^-a._ + au = A8 ^^-^^^^l;^^^^\ K'S). (264) 



-2(A; - A)a^^ + 4{k - 4)(Tm + (A; - 4)(A; - 5)aA,A + 2(A; - 4)aA,n 

0(2 
A;6 



+2an - 4^3 - 8au + 4a^ = 144 (^ 3)(2fe 3) ^ (^i,i,2,2). (265) 



2(A; - 4)(T^^ + 2(Tn + 6(A; - 4)aM + 8(Tii,ii + 16(7^ + 8(A; - 4)cta,ii 

+2{k -4){k- 5)aA,A + 8au - 4au = -432 ^^^,7^^ , 

~ k° 

(Dli'S) (266) 



2a^ + (A:-4)a^^-2au + 2a^ = 48 ^^ 3)(A: 3A: + 3) ^ 



^1112n 



(267) 



(A; - 3)(A; - 4)a^/ + 2(A; - 3)ct> + 4(A; - 3)au + 6a^ 
+24(Ti + 4ct« + 2{k - 3)0-^ + 2{k - 3)cx^ 

These coupled equations are equivalent to the following equations, 

,(2A;-3)2 



2aii,ii + (Tn = 36 

C^IIIJ + C7u = 16 



A;6 ' 
(A;2-3A; + 3)2 



-(/c - 3)o-in,i + 0-A 



A;6 

(A;-3)(A;2-3A; + 3)2 

^"^ k^ ■ 



a^+{k- 3)au = 8—{k - 3)(P - 3k + 3)^. 

rv 



(269) 
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From (268), we obtain by the use of (267), 

(A; - 3)a-ii + a a = 8^ ^-^— ^ (270) 

which is same as (269). 

We have the following equation by considering -D1233) 

6 (2k — 3) 

-(7A + an + 2a> + {k- 4)aM + 2a^ = 72^-—^ (271) 



k-3 -' ^ ' - ' ' ' ''' ' - k^ 

We find that many coefficients are not determined uniquely. We have 23 
coefficients and 10 coefficients are determined uniquely. 13 coefficients are 
not determined uniquely. There are 9 linear independent relations among 
these undetermined coefficients. These 9 linear independent relations are 
(259), (260), (263), (264), (262), (245) ,(249), (251), (255). We find there are 4 
quartic identities in appendix C, (282)~(285). Therefore, we are able to 
choose freely 4 coefficients. 

For instance, if we choose aM = — 36(A; — 3)^, we have 

au = -36(A; - 3)^ 

aA^n = l08{k-2){k-3) 

ail = -24(A; - 3){k'^ - 3k + 3) 

aA,i = -18(A;-l)(2A;-3) 

ai^^ = -12{k-3){P -3k + 3) (272) 

Appendix C: The cubic and quartic identities 

There is one cubic identity in the case k = 4 for the variables Tjj = 

(xj — Xj){\i — \j) as shown in [5], 

^3 = KV34] - [ri2ri3T24] + [^12^-137-23] = 0. (273) 

which may be conveniently depicted graphically by 

/3 = [II,I]-[N] + [A]. (274) 

For k > 4, this cubic identity is modified by a factor (A; — 3) as 

-^3 = [nV34] - [ri2ri3T24] + {k - 3) [ri2ri3r23] = 0. (275) 

We apply the differentiations D'^l^, where the indices take one of values 
among (l,2,3,...,k). By the symmetry, it is enough to consider the values 
(1,2,3,4). We find all possible such differentials satisfy this identity, 

D'aiTh = 0. (276) 
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This is a proof that we have an identity of (275). 

The geometric proof is also possible. In the case k=4, we consider the 
coincidence of the variable 0:4 -^ Xi and A4 — > Ai. In this limit, we obtain 
three type 3-point graphs. From [II, I], we obtain [Z.], and this [Z] is exactly 
cancelled by the same graphs generated from [N] in this coincidence limit. 
The graphs [N] also generate the triangle graphs [A], which are cancelled 
by the existing triangle graphs of the last term of the identity J3. Thus we 
prove that Is = for /c = 4. For A; > 4, we can prove this identity by the 
inductive method, namely for k = 5, when we take the limit X5, A5 -^ xi, Ai, 
the graphs reduces to the graphs of A; = 4 by the cancellations. 

Since we have this identity, the coefficients o"n,i, o"n, o"a are not uniquely 
determined as we have seen in the appendix B. 

In the fourth order (4 line graphs), the identity Is = at k=4 is general- 
ized by the multiplication of [T12] = T12 + T13 + ■ ■ ■ + r3_4. By writing all types 
of graphs, we find easily the following identity of k=4, 

h = hx [ru] 

= [III,I]+2[II,II] -[>] + [ A ]-4[n]-[ n ] = (277) 

This identity of /4 = at k=4 is verified directly by the differentiations 
D'^tT (a,b,c,d,k,l,m,n are 1,2,3 or 4). 

For the general k (A; > 4), we use the cubic identity of (275). Multipling 
[ti2\ on each terms of the cubic identity, we get (the indices are now from 1 
tok), 

[II, I][ri2] = 2[II, I, I] + [III, I] + 2[II, II] + [Z, I] + [3] + 2[A, II] (278) 

-[N][ri2] = -2[M] - [N,I] -2[>] - [H] - [D] - 4[n] - 2[ZZ] (279) 

(k - 3)[A][ri2] = (A; - 3)[A, I] + (k - 3)[ A ] + (k - 3)[>] (280) 

Adding these three equations, we obtain an identity. 

2[II, I, I] - [N, I] + (k - 3)[A, I] + [III, I] + 2[II, II] + [Z, I] - 2[M] 
+2[II,A] + (k-5)[>]-[ n ]-4[n]-2[ZZ] + (k-3)[ A ] = 

(281) 

This long identity is devided into 4 sub-identities. 

By the analogy of the cubic identity, we extract from (281) an identity, 

2[II,I,I]-[N,I] + (k-5)[A,I] = (282) 

This identity is the cubic identity plus a separate line. The factor (A; — 5) 
means that when A; = 5 all terms are vanishing. ([A, I] is nonvanishing. 
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therefore it should have a factor {k — 5) as a coefficient, since other two 
terms do not exist for k=5.) 

This simple (fundamental) identity at fourth order is proved by all possi- 
ble differentiations D^''" . For instance, with the differentials of DlH^, we find 
D\lll[II,I,I] = 4(k-3)'(k-4)(k-4), Dl\il[N,I] = 16(k-3)(k-4)(k-5) and 
D^iiiA,!] = 8(k-3)(k-4). These values satisfy (282). All differentiations 
which appeared in appendix B satisfy this identity. 

The second identity in (281) is a generalization of (277), obtained by 
adding a factor (A; — 3) for [ A ], 

[III, I] + 2[II, II] - [>] + (k - 3)[ A ] - 4[n] - [ n ] = (283) 

This identity is verified by the differentiations of D^l^J^, where the indices 
take values from one to four. Indeed in appendix B we have evaluated various 
differentiations, and we have found the contributions of every graph which 
enters in this second identity. We have checked all possible differentiations. 
As third identity in (281), we find 

2 [A, II] - 2[M] + 4[A, I] + 2(k - 4) [D] - (k - 4) [II, II] = (284) 

This identity is also proved by the consideration of all possible differentia- 
tions. 

Subtracting these three identities from (281), we obtain the fourth iden- 
tity, 

[Z,I]-2[ZZ]-2[A,I] + (k-4)[ > ]-2(k-4)[n] + (k-4)[II,II] = (285) 



We have checked this identity by the evaluation of all differentiations 

DS" in (247). 

Further, we check them by the differentiations -D^^^", (a,b,c,d,k,l,m,n = 

1,2,3,4), which appeard in the calculations in appendix B. They are Dllu^^luh 
P)U22 nll22 niii2 nll33 nll33 nll33 niiis niiis 

-^^1234 5-^^2234 )-'^2433 5 -'-^2233' -^^2244' -'-^2234'-'^2223'-'^2224- 

In the appendix B, we were left with 13 undetermined coefficients C[graph] 
. We have 9 linear independent equations for them here. Since 13 - 9 = 4, 
given that one has 4 quartic identities (282) ~ (285) all the coefficients may 
be determined at the expense of introducing 4 arbitrary parameters. 

In the case [3 = 4, there are no double or multiple line graphs. Therefore, 
we have no cubic or higher order identities, and all coefficients of the r ex- 
pansion are uniquely determined without ambiguities. In the next appendix 
D, we discuss the case of /3 = 4 in a more systematic way. 

Appendix D : Characterization of the r expansion by specific 
differentiations 
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As discussed in appendix B, the coefficients of the r expressions for the 
products of symmetric functions S4(x)s4(A) are determined by focusing on 
the monomials x^^x^^ ■ ■ ■ A^^AJ^ ■ ■ ■• Such terms may be selected with the help 
of the differentiations -D*^*^'.'.'.. Many terms are thereby uniquely determined 
since they are characterized by the existence of particular combinations of 
Xi and Xj. Since each r term has a specific topological structure when rep- 
resented graphically, one may characterize such r terms by the differentia- 
tions D*^*^'.'.'.. For instance, the term [rf2] = J2i<j'^ij is uniquely character- 
ized by xfX"^, which is equivalent to D^. Other types of r-terms, such as 
[ti2''"i3],[ti2T34], have no such x^Ag factor. Therefore xlXl is a unique factor 
which selects [T12]. 

Indeed we have used such properties for the determination of o", which is a 
coefficient of the r expression for S2{x)s2{X). We study such correspondences 
between the unique differentiations and the r terms (graphs), because we 
can apply these correspondences directly to determine the coefficients C for 
the HIZ integrals. 

This one-to-one correspondence between r terms and differentiations D 
may be used as follows. We write the function / in (2) as a sum, 

/ = 1 + /(I) + /(2) + /(3) + • ■ ■ (286) 

where f\i) is the 1-th order term of the r-expansion. The function / is given 
by 

which is expanded as 

/<.. - ^(i:.,)-+ ,(,4(,_,) ^3WMA) (288) 

Table D-1: Characteristic differentials of r terms and / 



1=1 [I 


D2[I] = -1, 




D^(f(i)) = i 


1=2 [a; 


D2MA] = 2, 




D23(f(2)) - k(k+a) 


II 


D^M = 4, 




Tillf'f A 2(1+q) 
J-'22lf(2)j - k(k+a) 


[i>i] 


Bll[l,l]=2, 




D34(f(2)) = k(k+a)(l 1 k-l) 


1=3 [I, I, I 


Dii[i,i,i] = 


—6, 


T-vl23/'f, .\ _ 6 ri 1 3a 1 2a'^ 


-L>'456lndjJ k(k+Q)(k+2Q) ^-^ ' k-l ' (k-l)(k-2) 


[A, I] 


r)ll2rA T] _ 


-6, 


^345(^(3)) k(k+Q)(k+2Q)(-'- ' k-l) 


[Y] 


■'-'234 [ J- . — 


-6, 


D234(f(3)) - k(k+a)(k+2a) 



When the r-terms (graphs) are selected by a single differentiation D, the 
coefficients of the corresponding term in the r expansion are given by 



D [graph] 



C[graph] — ^f ^_ ^i (289) 



For instance, from the above table, we obtain 

The values of the coefficients Cigraph] coincide with the values in Table B, 
when we put a = 2/(2 — /?). 

The third order term /(3) in (286) follows from the Jack polynomial ex- 
pansion, 

+ 3(A; + a)(fc + 2a)(A:-l)(A:-2) ^^(^^'^^(^^ ^^91) 

The coefficients of [II, I], [N] have arbitrariness due to the cubic identity, 
which has been discussed in the appendix C. We take here the reasonable 
constraint that the double line graphs have ^-^ as overall factor. This means 
that Cn,i is proportional to ^-^, and it is vanishing for a = — 1(/3 = 4). For 
the single line graphs, their coefficients are simply — p- in the large k limit, 
for the third order terms. With these constraints, we have from the table 
D-1, 

^"'^ = -2k{k + a){k + 2a)^^ + k^^ ^^^^^ 

^'' = ~k{k + a){k + 2a)^^^k^l^ ^^^^^ 

which coincide with the values in Table B by the substitution of a = ^^. 
By the differentiation -D33I, we have from Table D-1, 

-4Ca = -4(A; - 3)ai,i - 8C^^, + Dl^^.if^,)) (294) 

Using the value of Cn,i in (292), we obtain 

1 a^ 

^^ = ~k{k + a){k + 2a)^^~r^^ ^^^^^ 

which coincides with the value in Table B. 
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For the fourth order, we have from the expression of Table in the appendix 
A, 

" 3(A: + a){k + 2a){k - l){k - 2) % ^^^>3(£)33(A) 

a2(l8a2 + isak - ISa^A; + 6^ - ISaP + Ga^p _ 5k^ + SaA;^ + A;^) 
^ 8A;(A; + a)(A; + 2a)(A; + 3a)(A; + a - l){k - l){k - 2){k - 3) 

a^k{k'^ + ak — k + a) ~ 

^ 4{k + a){k + 2a){k + 3a)(A; + a - l){k - l){k - 2){k - 3)^4(a:)s4(A) 
a3(2A;2 + 3aA; - 3A; - 3a) 
4(fc + a){k + 2a){k + 3a){k + a- l){k - l){k - 2){k - 3) 
X [S2(£)%(A) + S4(5)(S2(A))'] (296) 

If we apply -D2345 '^'^ /(4)) *^ obtain 

24 
^^^^^f^^^^ = k{k + a){k + 2a){k + 3a) ^^^^^ 

This -D2345 is ^ characteristic differentiation of the r term [X], namely it is 
the unique differentiation for [X] term. Since we have -D2345[X] = 24, one 
obtains 

^"^ " k{k + a){k + 2a){k + 3a) ^^^^^ 

which coincides with the value given in Table B, if we put a = ^^. Other 
differentiations of r terms and /(4) are represented in the following table D-2. 
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Table D-2: Characteristic differentiations of r terms and /(4) 

1=4 [X] D^^|^[X] = 24, D^J^^(f(4)) = k(k+a){k+2a)(k+3a) 

[Y, I] D345g[Y,l] = 24, D345g(f(4)) = ^^Zr)(k+a)(k+2aKk+3a) 

W Al r)1122rA A] _ o/i nlll2/f \_ 24(fc+2a-l)(fc+3a-l) 

[il,iVJ Xy345g[YV,iVJ - Z^, -t^3456U{4)J - fe(fc+a)(fc+2a)(fc+3a)(fc+a-l)(fe-l) 

[A, 1,1] Dl^23[A,l,l]=24, 

r)U23( f \ _ 24 /i _ 2a , _7a_ , 6a^ 20^ 

-^4567U(4); fc(A:+a)(fc+2o)(fc+3a) V fe+a-l "*" fc-1 "•" (fc-l)(fc-2) (A;-2)(fc+a-l) • 

[1,1,1,1] Dii^[l,l,l,l] = 24, 

1234// \ _ 24 /I I 6q I 9a2 So^ 



7-)1234/f \ _ 24 (-, _|_ 6a , 9a^ , 

-^5678U(4); — fc(fc+a)(fc+2a)(fc+3o)V-^ "^ fc+a-1 ^ (fc-l)(fe+a-l) "^ 



fc(fc+o)(fc+2o)(fc+3o) V ' fc+a-1 ' (fc-l)(A;+a-l) ' (A;-2)(fc+a-l) 
I 25a3 8af , Ga" ^ 



(fc-l)(fc-3){A:+Q-l) (fe-2){A:-3)(fc+a-l) ' (fc-l)(fc-2)(fc-3)(fc+a-l) ' 

[^] ^2244P] = 48, ^2244(/(4)) = k(kTa)(k+2§)(kT3aj (^ + k^) 
[<] -D2233i<] = 96, -D2233(/(4)) = fc(fc+a){fe+2a){fe+3a) 
[H -D2234[H = 48, -D2234(/(4)) = fe(fc+a)(fc+2a)(fc+3a) 
[ ^ ] -^2223 i ^ J = 144, -D2223 (/(4) ) = fc(fc+a)(A;+2a)(fc+3a) 

[1111] D--[iiii] = 576, DMM(f(4)) = t:sii:sii:s 
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Table D-3: Characteristic differentiations of r terms and /(4) 

^33l5(/(4)) = fc(fc+a)(fc+L)(fc+3«) [(^ + ")(1 + |^) + (^ + ^)] 

^2ii(qA,II] [A, II] + C[M] [M] + C[A,A] [A, A]) 

= -24qA,n] - 2AC[M] - 12(A; - 5)qA,A] 

f)1122(f \ _ 12{k+2a-l){k'^+4ak-3k+3a'^-9a+2) 
-^^2345 U (4) j — k{k-l){k+a-l){k+a){k+2a){k+3a) 

^222l(C[iii,i][III,I] + C[u][n]) = 36C[in,i] + 36C[u] 

nlll3/'p \ _ 12{l+a){ak+2k+3a^+2a-2) 
-^2224l-^(4)J — k{k+a){k+2a){k+3a){k-l) 

/^2'2ll(C[ii,ii][II,II] + C[n][n]) = 32C[n,ii] + 16Cp] 

nll33/x N _ 8(l+a)2 /-, , J^ , 2a2(2+a) x 

-^2244U(4); — k(k+a){k+2a){k+3a)y^ ^ k~l "^ (l+a)(fc-l)(fc+a-l) / 

^ k{k+a)(k+2a)(k+3a)\^ ^ fe-l/ 



^222l(%]y + C[iii,i][III,I] + C[A][A]) = -72C^-36(k-3)Cni,i + 36CA 

nlll3/'/ \ _ 6(l+a)(2afc^+fc^+6a^A:-afc-5fc-12a^-7a+4) 
-^^22231/(4) J — fc(fe+a)(A:+2a)(A:+3a)(A:-l) 

/^ILKCpi,!,!] [II, I, I] + C[Y,I] [Y, I] + C[N,I] [N, I] + C[A,I,I] [A, I, I]) 

= -24q„,i,i] + 12qY,i] - 48qN,i] - 12(A; - 6)qA,i,i] 

nll23/'f \ _ _19 1+^ fi _ 2a _|_ 7a , 6a^ 2a'^ 1 

-^1456U(4); ^^(k+a)(k+2a)(k+3a)\-^ fc+a-1 """ fe-l "^ (fc-l)(fc-2) (fc-2)(fc+a-l) J 

48 M I 8a _ _4a of , (l+a)(4afc+3a^-4a) i 

fc(fc+a){fe+a)(fc+2a)(fc+3Q) ^ """ fc-1 fc-2 (k+a-l)(k-2) "•" (fc-l)(fc-2){A;+a-l) J 

^223l(q>][>] + c^iiMi] + c^m + q^lH) 

= -24q>] - 12(A; - 4)^/^] + 24Cy - 24C[^] 

nlll2/f \ _ 12(fc2+2afc-3A:-3a^-7a+2) 
-^^2234 U (4) J — k{k+a){k+2a){k+3a){k~l) 

DlllUquMl] + qM][M] + C[N,i][N,I] + C[A,i][A,I] + C[A,A][A, A] 
+%i][Z,I] + C[A,ii][A,II]) 
= 4q^^] - 24qM] - 8(A; - 5)qN,i] - 8qA,i] - 4(A; - 5)qA,A] - 8%i] - 8qA,ii] 

^2^it(/(4)) 
_ 4(12~78a+120a^-3ea^-6a'^-36k+151ak-120a'^k+7a^k+39k^-90ak'^+27a^k'^-18k^+17ak^+3k'^) 

~ k{k-l){k-2){k+a-l){k+a){k+2a){k+3a) 
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The coefficients in Table D-3 have ambiguities due to the quartic identities 
discussed in Appendix C. For the coefficients Cu ji, C'r//], we have 

(299) 
If we assume that Cr^ n has an overall factor (1 + a), which implies that it 
vanishes for /3 = 4, we uniquely determine these two coefficients as 

1 + a 3q; 

^[^'^1 = 2A;(A; + a)(A; + 2a)(A; + 3a)^^ + I^^ ^^°°^ 

and 

1 2q; 

^^^ " k{k + a)(A; + 2a)(A; + 3a) ^^ ^ I^^' ^^°^^ 

These values coincide with Table B, and for /3 = 4, they are consistent with 
Table A. 

Subtracting the value of — 12(A; — 5)Ca,a, which was obtained in Table 
D-2, from -D234|(/(4)), we obtain 

„ „ 1 + a /_, a; , ,^ 3a , 
Ca,ii + Cm = ^,,, , ,,, , , ,,, , , ^ (1 + r^ 7)(1 + 



2A;(A; + a)(A; + 2a)(A; + 3a) A; + a - 1 k~\ 



A;(A; + a)(A; + 2a)(A; + 3a) A; + a - 1 k-l' 

(302) 

Normalizing again the double bond with the factor (1 + a)/2 which vanishes 
for /3 = 4, we obtain 

r — 1 + a a 3a (onr/^ 

^^'" " 2k{kTaj{k + 2aj{kT3a)^ + A; + a-l^^ + :fc^^ ^^^^^ 



and 



1 a 2a 

^^ = A;(A; + a)(A; + 2a)(A; + 3a)^^ + k + a-1^^^ + I^^ ^^°^^ 

We write -D2224(/(4)) ^^ 

1113 , _ 6(l + a)(l + 2a) ^^ , 3a 



^2224 (/(4)) = ,/. , ^N/. I o^\fU I Q^a (^ 



A;(A; + a)(A; + 2a)(A; + 3a)' k - I' 
18(1 + a) a 

^ k{k + a){k + 2a){k + 3a) ^^ ^ ^^^ ^^°^^ 



Then we find that 



^ _ (l + a)(l + 2a) ^^ 3a 

^"^'^ " 6A;(A; + a)(A; + 2a)(A; + 3a)^^ + I^^ ^^"^^ 
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and 

^- " 2k{k + a){k + 2a){k + 3a) ^^ ^ I^^ ^^^^^ 

We have from Table D-3, 

{1 + aY Aa 2a\2 + a) 

^n,II = 7771 — ; 771 — , r. \fi — , o \ U + 1 7 + 



Ak{k + a){k + 2a){k + ?,ay k-1 {1 + a){k - l){k + a - 1)' 

(308) 

and 

1 2a 

^° " A;(A; + a)(A; + 2a)(A; + 3a)^^ ^ I^^ ^^°^^ 

For the differentiation -D2223' ^^ ^^^ obtained values of C^, Cni,i, then we 
find uniquely, 

r = 1±^ ri-^^^ r3im 

- 2A;(A; + a)(A; + 2a)(A; + 3«)^ k - l' ^ ' 

For -D14II, we subtract the two known expressions of Cy,i and Ca,i,i, and 
obtain the following two coefficients, 

1 + a 2a 7a Qa^ 

^n,i,i = TTTi — ; — T71 — , o \/i — , o \ i^ 



{2k + a){k + 2a){k + "iay k + a-l k-l {k-l){k-2) 

2a^ ^ , , 

(311) 



(A;-2)(A; + a-l) 



1 8a 4a 

^N,I — 771 — ; 771 — , o \l^ — I o \ i-*- 



A;(A; + a)(A; + 2a)(A; + 3a)' k-\ k-2 

a^ (1 + a)(4aA; + 3a^ -4a) 



(A; + a-l)(A;-2) (A; - 1)(A; - 2)(A; + a - 1) 

^. The factor [1 - ^^^_^ 



(312) 



We have chosen the overall factor of Cn,i,i as '^^. The factor [1 — ,^"_. + 



7q I 6a^ 2a^ 



in Cii^ij, is the same as for Ca,i,i. 



fc-l ' (fc-l){A:-2) (fc-2)(fc+a-l)J 

For Dim, ^^ obtain after the subtraction of the known three terms, 

r = I n Q^(^-^) a r3i3^ 

- A;(A; + a)(A; + 2a)(A; + 3a)^ k-l ' ^ ' 

This expression becomes m._i)2(l._2) fo^ ct = —1, and it is consistent with 

-^2235 

coefficients, one obtains Ca,i as 



the value given in Table A. For -D22355 using the previously determined 6 



^ 1 a(a-3) a2(4a-3) 



A;(A; + a)(A; + 2a)(A; + 3a)' A; + a - 1 {k-2){k + a-l) 
a2(3a2-2a + 3) . 



(A;-l)(A;-2)(A; + a-l)^ 
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(314) 



This expression reduces to fc(-fcii)2Cfc„2'i« ^'^^ a = —I, and this agrees with the 
result given in Table A. 

We have thus obtained explicitely all the coefficeints up to fourth order. 
They satisfy a remarkable property. When the parameter a vanishes 

becomes unity, and / in (287) is simply given by 

f = e~^^^<^^''. (315) 

Then the coefficients are given by 

C={-iy^ (316) 

where/ is the order of the r term, and C the corresponding coefficient ; g is 
a degeneracy factorfor the multiple lines of the graphs. For the double line, 
g = 1/2!, and for the triple line, g = 1/3! etc. 

The pole terms with denominators {k — 1) ,{k — 2) ,{k + a — 1),..., which 
come from -^-j^ in (287) disappear in the expression of the coefficients in the 

limit a = 0. This leads to the fact that the pole terms -j^, -^' k+a-i ^'"^ 
always appear with a factor proportional to a or to a power of a. Then, the 
coefficients C may be factorized as 

where g is a degeneracy factor, and for the single multiple line graph case, it 

is given by 

n^~ii(l+ma) 
g = j 318 

q is the number of multiple lines. When there are n mutiple lines, (qi,q2,---,(ln) 
in a graph, the degeneracy factor g becomes 

9 = [[ -1 (319) 

For instance, in the case of Cnji, the degeneracy factor g is — — . 

The correction of order a in (317) is order of | in the large k limit. 
Therefore, in the large k limit, all the coefficients of the r-expansion are 
given by (317). 
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The second remark is about the correction terms to (317). Common 
corrections appear when graphs have common structures. For instance, star 
graphs have no corrections. As star graphs, we have up to the fourth order, 
Cj, Ca,Cu, Cy,Ciu,Ci_, Cx,C«,C|=,C/,Ciiii. They are given exactly by the 
leading term of (317). 

For the graphs with one separate line, the correction terms are agaion 

all the same if the remaining part is a star graph. For instance, we find a 

2a 
common factor for C\^i and Cnj, which is (1 + -). For Cy,i,C/ j,Cni,i, 

rh 1 

the common is (1 H ). 

k — I 
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